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Current-voltage (/-V) characteristics of PVF , films grown from the solution and characterized by having (a + f) crystalline 
modification have been studied as a function of temperature and thickness in the M-P-M sandwich configuration. The /-V 
curves show two distinct regions of conduction, viz. the comparatively low-field ohmic region with a slope ~ | and the high 
field non-ohmic region with a slope 2. The high field non-ohmic behaviour has been attributed to the tunnelling of charge 
carriers into or via the traps, the tunnelling process being a thermally activated one. 


1 Introduction 

The discovery' of strong piezo-electricity in 
uniaxially or biaxially stretched and poled polyviny- 
lidene fluoride (PVF) films has drawn considerable 
interest in recent years particularly with regard to the 
mechanisms? responsible for this effect. PVF, films 
possess many crystalline structures* ~* such as a polar 
form ($-form), a non-polar form («-form), a mixture of 

ar and non-polar forms (f + «-form) and a 7-form®. 
, the B-form has been found to be mainly 
present inthe piezo-electric fi films of PVF;. Murayama 
1c Hashizumi® have established that in PVF, films 
taining both the B- and the a-forms, the piezo- 
the pyro-electricity increase with the 


t ““ ae mi a 


ia aa eee 
’ 


i pitty 


3 Results and Discussion 

In order to see the morphology of PVF, films grown 
from the solution, the IR spectra in the frequency 
range 530-480 cm~' were recorded at room 
temperature and were compared with the reported! 
Murayama’s spectra. Curves 1,2 and 3 in Fig. 1 
correspond to Murayama’s IR spectra for PVF, films 
whereas curve 4 corresponds to IR spectra of our films. 
Murayama’? had found that in PVF, films possessing 


the «-form, two absorption peaks occur at 520 and 485 


cm! (curve 1) and in those having the f-form, the 
peaks occur at 500 and 485 cm “ (curve 2). However, 
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in films possessing the («+ f)-form, three absorption 
peaks, namely, at 520, 495 and 485 cm “* (curve 3) are 
observed. As seen from curve 4 (present study) three 
absorption peaks at 515, 495 and 485 cm ~* have been 
observed. On comparing curves 3 and 4 (Fig. 1), it is 
seen that the films studied by us have absorption peaks 
very close to those reported by Murayama’ * for the (« 
+ B)-form indicating thereby that our films of PVF, 
consist mainly of the («+/)-form. Murayama and 
coworkers!?’!* have suggested that films having sucha 
crystalline modification would provide trapping sites 
and their electrical properties will be governed 
predominantly by the trapping mechanisms. Our TSD 
studies ’’* on such PVF, films confirm the dominating 
role of traps for the storage of charge; the charge 
carriers are trapped either at the microscopic or at the 
macroscopic level. The conductivity behaviour 
discussed ‘later in this paper also supports the vital role 
of the traps in governing the electrical properties of 
solution-grown PVF, films. 

Fig. 2 shows the room temperature /-V 
characteristics in log-log scale for PVF, films of 
different thicknesses. Curves 1 to 5 correspond to 
thicknesses 18, 27, 35, 57 and 80 wm, respectively. It is 
seen from these /-V curves that there are distinctly two 
regions of conduction, viz. the comparatively low-field 
ohmic region with slope ~1 and the high-field non- 
ohmic region with slope <2. Further, the transition 
voltage (V,,..) at which the non-ohmic behaviour starts 
appearing is seen to depend on the thickness of the 
films; it increases with the increase of thickness as is 
-* in the inset of Fig. 2. 
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Fig. 2—/-V characteristics of pure PVF, films a 

(303 K) for different thicknesses (Curves | to § co 

thickness of 18, 27, 35, 57 and 80 am, respectively nine : 
shown a log-log plot of Vern V8 d] ne 


{room temperature 
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Fig. 3—/-V characteristics of pure PVF, films (thickness ~ 35 um) 
for different temperatures [Curves 1 to 7 correspond to 
temperatures of 303, 317, 323, 333, 353, 373 and 393 K, respectively] 


Fig. 3 shows the /-V curves in log-log scale for PVF, 
films (thickness ~35 ym) at different temperatures. 
Curves | to 7 correspond to the temperatures 303, 317, 
323, 333, 353, 373 and 387K, respectively. It is seen 
that at low temperatures the same two distinct regions 
of conduction as stated above are observed whereas at 
high temperatures, there is practically no break in the 
I-V curves. Thus with the increase of temperature, i.e. 
as One goes from curves | to 7, the slope of the high- 
field conduction region decreases. 

The high-field electronic conduction mechanisms 
commonly discussed for various polymer films are 
space-charge-limited conduction (SCLC), Richardson- 
Schottky (RS), Poole-Frenkel (PF) emission and 
tunnelling of carriers into or via traps'*"5. In order to 
establish the dominant mode of conduction in a 
particular material, one has to look into the detailed 
analysis of /-V data of the material in terms of the 
theoretical considerations available for different types 
of processes. : 

For the currents to be space-charge-limited in the 
high-field region, the current density J in case of 


shallow traps should obey'® rere eel 
2 HeOV? y _ 8enod? nap ean i ai 
8d “an 980 
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where yw is the drift mobility, ¢ the high frequency 
dielectric constant, V the applied voltage, V,,,, the 
transition voltage separating the two regions of 
conduction, d the thickness of the film, e the electronic 
charge and @ a constant which describes the effect of 
traps and is the ratio of the free carrier density (ng) to 
the trapped carrier density (n), i.e. 9=no/n.,. 

It is seen from Eq. (1) that (i) log-log plot of / vs din 
the V? region should be a straight line with slope —3, 
and (ii) log-log plot of V,,,, vs dshould bea straight line 
with slope 2. Such plots were made and the slopes of 
—2.2 and 1.2, respectively were obtained which rules 
out the possibility of space charge limited currents. 

If one considers the RS or PF emission as the 
dominant process of conduction then the current 
density / follows the relation'*’'* 


T=Iyexp(BE'?/kT) ber tt) 
where J)>=AT*exp(—/kT) for RS 
and J)=0,V/d for PF 


and the constant 
Brs =(e?/4 TEE)! !? = Bor 


where &) is the permittivity of free space, A the 
Richardson’s constant, do the low field conductivity 
and E the applied field, ic. E=V/d. Eq. (2) suggests 
that the plot of log/ vs V'/? should be a linear one 
having a slope of B/kT. Such a plot was made for PVF, 
and is shown in Fig. 4. The value of B was estimated 
from the straight portions of the graphs and the value 
of high frequency dielectric constant was calculated by 
using Eq. (3). It was found that on substituting the 
value of Bas, €= 13 and on substituting the value of Bpr, 
e252. However, as the reported’” value of ¢ for PVF, 
at 10° Hz is ~2.98; such a large discrepancy in the 
estimated and the reported values of € completely 
eliminates the possibility of RS or PF mechanisms 
r the high field currents. 
Now, the /-V data on PVF; is analysed i in terms of 
ee seeekine of carriers into or via traps, i.e. traps- 
d tunnelling (TAT); the process, however, may 
be a ‘thermally activated one. The Danese 
expe Bod current density 7 in the case of 
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Fig. 4—1-V data of Fig. 2 plotted as log/ vs V'/? [Curves | to 5 
correspond to thickness of 18, 27, 35, 57 and 80 um, respectively] 
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Fig 5—Plot of activation energy (U) vs the sq. root of applied voltage 


determined from the conventional Arrhenius plot of 


, log/vs 1/T will bea field-dependent one. This indeed is 


ee A plot of activation energy vs sq. 
the applied voltage is shown in Fig. 5. It is seen 


4 os Tao ‘from 1.02eV at 3V to 0.72eV at 


ees 


pee 


| fal “325 
. . 
os “—. o . 


INDIAN J PURE & APPL PHYS, VOL 21, JUNE 1983 


1000 V. Such a field dependence will especially be 
noted in the TSD studies’. Our TSD studies’ on 
PVF, support this, as the activation energy of the 
charge carriers responsible for the TSD relaxations has 
been observed to depend on the polarization field. 
Therefore, at high fields, the trap-assisted tunnelling 
into or via traps may be suggested as the dominant 
mechanism of conduction in PVF ,; the process, 
however, being the activated one. 

The effect of temperature on the /-V characteristics 
may now be explained as follows. With the increase of 
temperature, the probability of thermal ionization of 
the trapping centres increases thus causing a shift in the 
quasi Fermi level which gives rise to a lowering of the 
barrier across which the electrons have to be 
transported and the conduction becomes more or less 
ohmic. Further experiments are in progress to study 
the effects of dopants and electrode material on the 
conduction behaviour of these films. 
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Thin-film components are useful for miniaturization and high frequency performance of electronic circuitry. The present 
paper describes the high frequency performance of thin-film capacitors using composite films of SiO + 5% B,O3. These films 
yield thin-film capacitors having low losses and are useful for high frequency applications. The performance of these capacitors 
is also compared with that of SiO thin-film capacitors. It has been observed that at high frequencies the electrode parameters 


also contribute towards the losses. 


1 Introduction 

SiO is commonly used as a dielectric for thin-film 
capacitors, though the iosses of such capacitors 
increase at high frequency. The dielectric loss of this 
material depends on the processing parameters! ° as 
well as on doping of the material*” ©. The addition of 
impurities like B,O, [Ref. 5,6] to SiO reduces the 
dielectric losses to a great extent due to increase in non- 
crystallinity of the film’. But it has been observed that 
the electrode parameters affect the dielectric losses of 
the films at high frequencies*’’. In the present paper, 
the variation of capacity and dielectric losses of thin- 
film capacitors of composite film of SiO+5% B,O, 
have been compared with the variation of these 
parameters with thin-film capacitors of SiO film only. 
The frequency range of interest has been restricted 
from 1 to 100 MHz. The behaviour at lower frequency 
j is not of much interest as no large difference in 


dielectric properties of the materials is expected. 


3 2 ne etna Details 
are fabricated by sandwiching the 
films between the two conducting aluminium 
The electrodes are deposited on glass 
page evaporation of aluminium 
mask. Dielectric material is 
are nbeller'® source heater through 
tes types of dielectric materials have 

ite i. aad ots (SiO +5% B,O3). The 
n kept at 250°C. Three 


All the capacitors made in a set have a common 
electrode for ease of measurement; such a set of 
capacitors is schematically shown in Fig. 1. 

Electrodes are deposited simultaneously while 
fabricating capacitors with two types of dielectric 
materials, so as to ensure the same thickness and 
conductivity of the electrodes. This makes it possible to 
compare the behaviour of two types of dielectric 
materials. The thickness of the dielectric material is 
kept at 5000 and 7000 A as measured with the help of a 
quartz crystal monitor (model DTM-3, Sloan, USA). 
Other conditions of evaporation like substrate 
temperature, position of substrate with respect to 
quartz crystal monitor, source temperature etc. of 
dielectric materials, are kept the same for both 
materials. 

The capacitance values are measured for the 
capacitors prepared in the configuration of Fig. 1, 
using a Wayne-Kerr Bridge, model B-80 (V.H.F. 
admittance bridge). The capacitors have area ranging 
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from 0.5 to 4mm?. The capacitance of the smallest one 
is about 50 pF for 5000 A thickness of the dielectric. 


3 Theoretical Model of Capacitor 

The theoretical model for the estimation of capacity 
and dielectric losses is presented first before discussing 
the experimental observations. The model is based on 
the analysis done by Prozorvoskii and Negodenko'’ 
for dielectric losses only. 

Consider a thin-film capacitor (Fig.2) having a 
dielectric of thickness d. Assuming that there are no free 
charges between the electrodes, the potential field will 
satisfy the Laplace equation. Let d be much smaller as 
compared to the length / and the breadth b of the 
capacitor. Then the field ¢ will satisfy the two- 
dimensional Laplace equation: 

Cp 
Ht geno .. (A) 

If a potential difference V is applied across the 
electrodes, having conductivity a, at points having 
coordinates (1,+d/2) and (1, —d/2), the solution of the 
equation is 

_V sinay cos hax 

~ 2sin (ad/2)' cos hal 
where « is an arbitrary constant (to be defined later to 
satisfy the present model). The surface charge density 
on the metal film (t) is related to current density by 


do 
aoe aod | ee to 


ertte) 


where @ is the angular frequency. 
The current density in the plate is 


a OY 
ae . (4) 


y 


7 


i | 
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Fig. 2—Model of a capacitor 
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The distribution of charge at the boundaries of the 
dielectric (y=d/2) is associated with the normal 
component of the field density by the expression 


‘ 
t= 4. gh y= +d/2 =e (5) 
oy 
where € is the dielectric constant of the film. The 
positive value of charge refers to the upper plate and 
negative value to the lower plate. 

Eqs (3)-(5) yield, 


42 
+= pie ... (6) 
Ox oy 
From Eggs (2) and (6), 
ga en igl al 
o 2 
If aS <I, cot aS ~2/ad 
Therefore, 
ime 1 
~eaPD 
or 
2iwe 
=4+ /<— 
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The presence of free charges in the dielectric is easily 
taken into account by assuming the dielectric constant 
to have the complex value 

= &' —jg/w —je" ee 
where é’is the real part of the dielectric constant, g the 
volume conductivity of the dielectric and &” the 
imaginary part of the dielectric constant, which 
characterizes the losses due to relaxation. 

For SiO or a mixture of SiO and B,Os, the electron- 
relaxation and_ ion-relaxation polarization are 
characteristics of the material and, in practice, the 
value of &” does not depend on frequency in the rf 
range!?. Substituting Eq. (9) into Eq.(8) and taking 
into account that in the usual frequency range for a 
Capacitor, €”/e’ <1 and g/we’ <1, one obtains 


g g 


— e” +2. oY 
= , : vT' WE @ . 
% = a + jax - SI ( tae +i(1 er 
+ (10) 
The current through the capacitor, 
V 
l=i =e ; 
i()b 5 aob tan hal «. (11) 


Here i(/) is the current denne’ at a point ha 

coordinates (I,d/2). The negative sign aie right side 
of wh indicates that the current flows in the 
Opposite direction of x-axis. It has ‘ignored it . 
subsequent calculations. me seri ger 
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The complex admittance of the capacitor is given by 


I «aab 
Fa5= > tan hal usta) 


The dissipation factor, tandé, of the capacitor is 
determined by the ratio of the real part of the 
admittance to the imaginary part, 


Vreai 


tan d= . (13) 


By substituting Eq. (10) into Eq. (12), separating the 
imaginary and real parts, expanding the trigonometric 


and hyperbolic functions in Maclaurin’s series, and 


neglecting higher order terms, one obtains, 
: é” ® 21? we 

y-- Sed 
The first term in this equation is characterized by 


losses due to conductance. It is large at low frequencies. 
At high frequencies, the third term, governed by the 


tand= 5 .. (14) 


losses in the plates, predominates. 


Consequently, in capacitors having large ! values, 


‘the losses due to electrode parameters will be 


dominating, while in the smaller capacitors the losses 


_ due to dielectric relaxation will dominate. 


~ Se sgigagil is given a 


eS i] “és Bers es a i 
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noted that because of the factor |* the values of the 
dielectric losses and capacitance are dominantly 
influenced by the second term involving /?. 

For large area capacitors the effect of losses due to 
conductance and relaxation are almost completely 
obscured by losses due to electrode dimensions. 
Consequently, the measurements have been restricted 
to the smallest area capacitors (0.5 mm?) alone so that 
the dominance of the term due to electrode parameters 
in the Eqs. (14) and (19) is avoided (to estimate the 
conductance and relaxation losses). The effect of 
thickness of the dielectric (d), which only moderately 
affects the magnitude of this term has, of course, been 
ignored. In fact, since it occurs in the denominator of 
the term, an increase in its value further tends to offset 
the dominance of /?, occurring in the numerator. 

Fig. 3 shows the variation of capacitance with 
frequency for SiO of thickness 5000A and (Si0+5% 
B,03) for two different thicknesses, viz. 5000 and 7000 
A respectively, while the curve III is the theoretical 
curve for the same composition and thickness 7000 A, 
computed by using Eq.(19). The curve IV is the 
experimental curve for SiO films of thickness 5000 A.A 


comparison of the curves I and 1V,which correspond to 


the same film thickness but different materials, shows — 


that the variation of capacity with frequency for SiO is 


faster than for (SiO +5 % B,O3). The value of €’ ", which Bese ns 


represents losses due to relaxation of the diel 
ial, is not exactly known, as most of fie 
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making use of the same value in Eq. (19), the theoretical 
curve 1] has been plotted for a film thickness of 7000 A. 
This curve has a close fit with the experimental curve II 
obtained for the same film thickness. It follows that the 
value of €” computed from curve I represents the 
correct value. Further, since the slope of curve I is less 
than that of curve IV, the value of &” for the films of 
(SiO +5 % B,O;) turns out to be smaller than for SiO. 
It implies that use of the former as dielectric material 
should improve the performance of the capacitor at 
high frequencies. At higher frequencies, the inductance 
offered by the thin-film leads of the electrodes of the 
capacitor deserves special consideration. The 
capacitance is measured with an admittance bridge 
and at high frequencies the impedance offered by the 
inductance may be large enough of offset the effect of 
impedance of the capacitor, as follows. The value of 
inductance'* per unit length is given by 

b+h 


3 1193402245") 


i l 
L=5.08 x 10 Hl loses; ] 
non sh20) 


where b, h and ! denote respectively, the width, the 
thickness and the length of the electrodes in mils. If h 
<b, as is usually the case, the expression becomes 


b 
L=5.08 x 1071 loge} + 1193-+0.2264 | 


aes, (24) 


For typical values of /and b,as those employed in the 
present study, ZL has values around 1.0 nH/mm; 
consequently at high frequencies it may offer sufficient 
impedance to offset the effect of capacitance. The 
curves in Fig. 3 show that the capacitance decreases 
with frequency, at high frequencies of 10°Hz and 
above. From Eq. (10), it is seen that this results due to 
the effect of electrode parameters | and o. Based on the 
measurement of capacitance values, Hirose and 
Wada" have observed a decrease in dielectric constant 
at high frequencies around 1 MHz, which they have 
attributed to the dispersion. However, it is possible 
that the observed decrease results from the effect of 
electrode parameters as in the present work. 

Since the capacitance shows a wide range of 
variation with frequency as in Fig. 3 it is desirable to 
specify the values of frequency where the capacitance is 
reduced to approximately 3/4 of the Original value. 
These values lie at 20 and 60 MHz for SiO and SiO 
+5% BO; respectively, for the same thickness (5000 
A) of the dielectric material. Thus the addition of § yA 
BO, to SiO has upgraded the Performance of the 
capacitor to much high frequencies. 

The variation of tand at high frequencies j 
range 10-100 MHz has been shown in Fig 4 for the ie 
of (SiO + 5% BO) of thickness 5000 and 7000 A and 
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Fig. 4— Variation of dielectric losses with frequency 


for SiO films. A comparison of the curves shows that 
the addition of 5% BO; to SiO considerably reduces 
the dielectric losses. The addition of BO; leads to a 
suppression of crystallization of SiO and hence to the 
reduction of dielectric losses. The comparison of curves 
II and III for (SiO + 5 % B,O;) further shows that the 
dielectric losses decrease when the film thickness is 
increased. This is in agreement with the observation of 
Prozorvoskii and Negodenko!® who have reported a 
reduction in dielectric losses with thickness for SiO 
films. Later, Ravez'* has found that the increase in the 
film thickness increases the non-crystallinity of the 
amorphous films of SiO,. Therefore, it may be 
concluded that this should lead to a decrease in the 
dielectric losses as observed in the present work. 

Measurement of capacitance over a period of one 
year is made for a set of capacitors. These capacitors 
are found to be stable; the variation in the capacitance 
of a typical capacitor over a period of one year gives 
almost a straight line parallel to the time axis. This 
gives a stable value of capacitor. 


5 Conclusions 

The observed low dielectric losses at high 
frequencies for the (SiO +5 % B2O3) films suggest that 
capacitors with higher stability and suitable for high 
frequency performance can be made using (SiO +5% 
B,Os) as the dielectric material. . 

These losses can be reduced further if the electrode 
material offers low conduction losses, ie. it possesses 
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ee 
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higher conductivity. At high frequencies, the inductive 
effect of the electrodes also plays an important role. 
Therefore, while designing the thin-film capacitors for 
high frequency applications the electrode leads and 
_ plates should be smaller in length. 
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Normal coordinate analysis is made for the out-of-plane vibrations of ortho, meta and para iodonitrobenzenes using a 14- 
parameter modified valence force field. The force constants are refined by the method of least-squares employing the 
vibrational frequencies of the three molecules. The potential energy distributions are computed and the nature of the 


absorption bands is discussed. 


1 Introduction 

The vibrational spectra of mono substituted 

nitrobenzenes have been reported by Green and 

Harrison! with tentative assignments. Varsanyi* has 

summarized them along with the vibrational spectra of 

seven hundred benzene derivatives. The normal 

coordinate analysis of these compounds, which is very 

essential for a complete understanding of vibrational 

modes’, is yet to appear in the literature. In Part lV of 

this series, we have reported a modified valence force 

field for the in-plane vibrations of mono iodonitroben- 

a zenes*. This paper presents the results of normal 

coordinate Sait of Secor vibrations of the 
™ an ie 7 ms | 


1ate Treatment 

puting procedure is the same as that 
i our aise paper*. The structure and 
t-of-plane eee elie of the 


calculations converged (Only 30 fr 


parameters are those used earlier*. The orthonormal ‘ 
set of symmetry coordinates used in the calculations 
are given in Table |. The vibrational frequencies of o0-, 
m- and p-iodonitrobenzenes are taken from those — 
reported by Green and Harrison'. The initial set of 
force constants is obtained by transferring the ag 
potential constants from those of mono bromonitro 
enzenes°. In order to avoid large corrections at any 
time, initial refinements are carried out ne ade 
factor of 0.10. Only eight of the f “r 
constants are allowed to be refined ares 


compounds could be used for 
pate value for & <a torsi 


as ee ee ee 
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process). Two more interaction constants are fixed at 
values calculated with damping. The damping 
factor is then made zero and further refinement is 
continued. The calculations converged in two cycles 
under these conditions with an average error of 
93cm ~' between the observed and calculated 
frequencies. 


3 Results and Discussion 

The initial and final sets of force constants along 
with their dispersions are presented in Table 2. The 
observed and calculated frequencies and potential 
energy distributions (PED) in internal force constants 
for the three molecules under investigation are given in 
Table 3. The characteristic feature of PED in internal 
force constants for out-of-plane vibrations is that there 
will be large contributions, both positive and negative, 
from off-diagonal force constants apart from the fact 
that the contributions from certain diagonal constants 
exceeds 100% by large amounts®’’. As such, it becomes 
difficult to assess the actual contribution from a 
diagonal constant. This presents difficulties while 
interpreting the results. To circumvent these 
difficulties we have calculated the PED in principal 
internal coordinates by suitable addition of 
contributions from off-diagonal constants to the 
contributions from diagonal constants (A diagonal 
internal force constant which represents a group of 
equivalent internal coordinates is known as the 


Table 2—Out-of-Plane Force Constants (in mdyne A/rad?) 
of Mono Iodonitrobenzenes 


Force Coordi- Common Initial Final _ Dispersion 
Constant nates atoms value value 
symbol involved 
Yr =H _ 0.4530 = 0.4531 0.0009 
Yri er | — 0.6242 0.6275 0.0101 
é Cc-C = 0.3955 0.3959 0.0018 
Ym C-N _— 0.4059 =: 0.3917 0.0068 
Yai N=O — 0.2621 0.2629 0.0026 
C-—N _— 0.0001 0.0001 0.0 
C-H, 
C-H C-C -—0,0648 —0.0658 0.0 
c-C, 
Cc-C ccc 0.0103 0.0103 0.0 
¢ Ls iF e 
C-H C-C -—0.1085 —0.1085 0.0 
[a os N, 5 
C-H C-C -0.0311 —0.0311 0.0 
C= 
CoC C-C 0.2260 0.2260 00 
“4 @-C.. 01744... 01741... 00 


0.1878 


principal internal coordinate. For example, vg is the 
principal internal coordinate which stands for the four 
CH out-of-plane bends). These PEDs are also shown in 
Table 3. PEDs below 10% are not shown. The same 
table contains the vibrational assignments in Wilson’s 
notation’. 


3.1 Vibrational Assignments 

From Tables 2.and 3 it is clear that: 

(1) Observed and calculated frequencies agree well 
within the limits of experimental error and 
anharmonicity corrections, 

(ii) The final values of the force constants are not 
much different from the initial values as the 
corrections have never been very large, 

(ii) The dispersions of the force constants are quite 
small and 

(iv) The constants are comparable with those of 
related molecules. 


These observations indicate that the force field 
obtained in the present case is a physically reasonable 
set. Therefore, the assignment of various fundamental 
modes can be done on the basis of this vibrational 
analysis, rather than from hitherto qualitative 
considerations'’?. 

C—AH out-of-plane bending vibrations—In disubsti- 
tuted benzenes of the present type, there are four C—H 
out-of-plane bending modes designated as 5, 11, 17a | 
and 17b in the ortho and meta compounds and 5, 10a, 
17a and 17b in the para isomer. The corresponding 
phase relations in the eigen vector are: 

+=+-;+++4; ++-- and +--+ 
As identified from these phase relations and the PED 
presented in Table 3, the bands at 984, 777, 852 and 952 
cm —' in the ortho compound are due to the C — H out- 
of-plane bending modes 5, 11, 17a and 17b 
respectively. The corresponding vibrations in the meta 
compound are at 980, 798, 925 and 893 cm “' whereas 
these are identified near 965, 814, 841 and 954cm * in 
the para isomer. Thus, mode 5 with the phase relations 
+—+-— is the highest out-of-plane C—H bending 
vibration in the three molecules studied. This result is 
against the expectations of Green and Harrison! for 
para compounds wherein they infer 17a to be the 
highest C —H out-of-plane bending mode. The lowest 
out-of-plane C —H bending mode is that in which the 
corresponding internal coordinates either increase or 
decrease simultaneously. It is found that 17b is greater 
than 17a in ortho and para compounds whereas the 
reverse is true in the case of the meta isomer. This result 
agrees with the expectations of Green and Harrison’. 
Further, the present calculations show all the C—H 
out-of-plane bending modes to be pure as can be seen 
from Table 3. 
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Table 3—Potential Energy Distribution (PED) of Out-of- 


Plane Vibrations of /™ono lodonitrobenzenes 


Obs. Cal. PED in internal force constants PED in principal Vibrational 
freq. freq. internal coordinates assignment 
(cm~') (cm™') 
ortho \odonitrobenzene 
(a species) 
139 128 1071, 270, 287m, 227m, *° Yer. 4 18yp1, 450, 397m 10a + 16b + 10b 
194 207 10y_. 25yR1, 265, 627m, 307e1. é 40) pi, 216, 4970 10a + 16b + 10b 
410 416 —- 29yg, 11 ¥e1» 2159, 6470 1026 16a 
464 477 —- 12g, 87yR1, 1496 726, 34yR1 16b + 10a 
693 662 257g, 18yR1, 1798, 137 a1 11065 4 
728 734 10, 575, 297m, 737 a1 T3y m4, 186 o(NO,)+4 
Ti? 769 120yp 92yR 11 . 
852 863 1217p 1017, 17a 
952 950 123yp, 136, lly, Yr 109), 17b 
984 982 138yp, 300, 3272, Yr 120yz 5 
meta lodonitrobenzene 
(a” species) 
129 143 26ye1, 647R1, L1ypi, 5; 106 32yp1, 647m, 156 10a + 16b + 10b 
186 199 13y7p, 18yp, 416, 19yu, 297 R1> é, 107m.5 3371, 496, 247m 10a + 16b+ 10b 
414 416 = 3 lyg, 2160, 6974 946 16a 
466 471 21ypz, 66yR;, 1526 816, 30yR; 16b+ 10a 
663 . 669 14yp, 32ypR,, 1796 1146 4 =i 
728 728-255, 237m, 80Ye«1 807 m1, 126 o(NO,)+4 
798 791 1217p, 175 89), 11 x “ee 
p+ S99 886 =:126y, 126 1077p 17b 
925 927 121 yp, 226, 16yp;, YR 104) 17a 7s 
980 978 134yp, 265, 27yp, Yr 118), 5 i 
para \odonitrobenzene oo: aa 
(b, species ) a >: Lai 
96 + a : 
260 £ ‘ ey Brae page 227R1; 6, 18)m.s 2Sypri, 406, 387 . 10b+ 16b4+ 11 +% ry ake 3 " 
ee rey HE ee 371» 106, 567m 10b+16b+11 ce “=3 
ied - Sy, 67 ya, 1896, 527m 796 33y 1 _ “1 ‘a 
- 673 6582 : . » 9°YR1 6b + 10b vdiee e a 
: 267 x» 297p1, 1646, 14y m4, 996, 147 x44 = 
734 734535, 28 ym, 757m, 159.238 
65 SO ame ism ieee oan 
> oe hn Te: 13 187R, YR 1197p 
- 4 a ph . . ’ (a, species) 
—- og 440-33 yg, 1755 | pit 1208 
a wt 826-1189, 148, ria Yr * 
| — 838 ANTy iT 1007. nee 


iy ' a= ee a 


—- 5 - e — 


| 7 BRE bas Tae 


ee See To A ee 


eT ict 


ae dee ree enteeg —— acca we 
aap seers nthe re cies retoiag 
a a a 


RAO & RAO: NORMAL COORDINATE ANALYSIS OF IODONITROBENZENES 


of 16b in ortho and meta isomers whereas the 
assignment could be reversed in the para compound. 
The present calculations establish that frequency of 
16b is greater than that of 16a in the three isomers 
considered. Further 16a is a pure mode whereas 16b 
has a considerable mixing with C—I out-of-plane 
bending mode in the three molecules. 


NO, out-of-plane vibrations—Wagging and _tor- 
sional vibrations of the nitro group fall under this 
category. According to Green and Harrison*, @(NO,) 
is expected near 777, 728 and 734cm “' ino-, m- and p- 
compounds respectively. But, as explained earlier, the 
band at 777 cm~' is a C—H out-of-plane bending 
mode I! in the ortho compound. According to the 
present calculations w(NO,) is identified at 728, 728 
and 734 cm™ in ortho, meta and para iodonitroben- 
zenes. It is important to note that this mixes with mode 
4 in the three molecules. t(NO,) could not be included 


in the calculations as the corresponding experimental — 


peeaet available. rage it is not diseussed 


: a ers ies 


Sees . s eS 28 


a ne 
re ger: Fike sah vee 0 


(ii) The out-of-plane bending modes of C —I and C 
—N mix with each other with additional contributions 
from the ring torsion 16b. 

(iii) The order of the ring torsions is 4>16b> 16a. 
Further, modes 4 and 14a are pure whereas 16b 
invariably mixes with C—I out-of-plane bending 
vibration in the three molecules. 

(iv) Out-of-plane bending vibration of the nitro 
group mixes with ring torsional mode 4 to a small 
extent. 
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Normal coordinate analysis of ortho, meta and para bromonitrobenzenes is made using a 14-parameter modified valence 
force field for the out-of-plane vibrations. The force constants are refined by the method of least-squares employing the 
vibrational frequencies of the three molecules. The potential energy distributions are computed and the nature of the 


absorption bands is discussed. 


1 Introduction 

It is well known that, for complete understanding of 
the vibrational modes of fairly complex molecules, 
normal coordinate analysis is indispensable! *. For 
mono bromonitrobenzenes, even though the vib- 
rational spectra are reported with tentative 
assignments*°, their quantitative treatment using 
normal coordinate analysis is yet to appear in the 
literature. Hence, we made a rigorous mathematical 
study of fundamental vibrations of mono bromonitro- 
benzenes and report the results obtained on their out- 
of-plane vibrations in this paper. 


2 Normal Coordinate Treatment 

The computing procedure and the structure 
parameters are the same as those used in an earlier 
paper® with the C—Br bond length taken as 1.85A. 
The structure, definition of out-of-plane internal 
coordinates and the Symmetry coordinates are the 
Same as those used in part V of this series’. 


Table 1—Out-of- 


The frequencies of the out-of-plane fundamental 
modes of o-, m- and p-bromonitrobenzenes are taken 
from those reported by Green and Harrison*. The 
starting set of force constants is transferred from out- 
of-plane calculations for dihalobenzenes® (for the 
bromobenzene part) and nitrobenzenes? (for the C 
— NO, group). The C—N torsion constant is taken to 
be very low (0.0001 mdyne A/rad?) and fixed during 
refinement as the corresponding torsional frequencies 
are not available for the title compounds. The 
interaction constants for which no values are available 
are taken as zero initially. In order to avoid large 
corrections at any time, initial refinements are carried 
out with a damping factor of 0.10. Only 8 of the 14 
constants are refined in this way using 30 frequencies of 
the title compounds until the calculations converged. 
At this stage, two more interaction constants are fixed 
at their values obtained with damping, damping factor 
is made zero and the refinement is continued. The 
calculations converged in two cycles under these 


Plane Force Constants (in mdyne A/rad?) of Bromonitrobenzenes 
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Force Coordinates Common 
constant involved atoms 
symbol 
YR C-H — 
Yri C-—Br — 
C-C — 
Yu C-N — 
Yun N=0 — 
t C=-N is 
Yrs Yr C-H,C-H C-C 
95 Oi4, C-¢,¢-c Cc 
Pris Vr C-—Br,C-H C-C 
Yur Yr C-N,C-H C=¢€ 
Yr 6 C-Br.C-c C-¢ 
Yr 6 C-H,C-c C=C 
Yu. 6 C-N,C=c C-C 
Pus Yai C=N, C=Br C-¢ 


-—-—hUlUlhT TT 


Initial Final Dispersion 
value value 
0.4337 0.4530 0.0009 
0.5827 0.6242 0.0096 
0.3746 0.3955 0.0018 
0.5408 0.4059 0.0072 
0.6117 0.2621 0.0027 
0.0001 0.0001 0.0 
—0.0733 — 0.0648 0.0 
0.0103 0.0103 0.0 
—0.1085 —0.1085 0.0 
—0.0311 —0.0311 0.0 
0.2260 0.2260 0.0 
0.1741 0.1741 0.0 
0.0 0.1875 0.0 
0.0 —0,0463 
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Table 2—Potential Energy Distribution (PED) of Out-of- 
Plane Vibrations of Mono Bromonitrobenzenes 


Obs Cale. PED in principal inter- Vibrational 
freq freq. nal coordinates assignment 
coa™* so cm™' 
ortho Bromonitrobenzene 
(a species) 
132 131) 14 Yay, 45 4, 40 yg 10a + 16b+ 10b 
199 217 30 yp,, 44 Vag 10a + 10b 
413 419 1056 l6a 
462 479 716,37 yp, 16b + 10a 
690 661 111 6, 13 yay, 4 
730 735) 72 yu}, 176 w(NO,)+4 
175. -772—s«*OW Dp “1 
852 864 106 yp 17a 
953 949 109 pp 17b 
985 980 120 yp 5 
meta Bromonitrobenzene 
(a” species) 
147 148 64 yy, 21 yp; 10b + 10a 
193 206 34 yp;, 50 5, 13 yy 10a + 16b + 10b 
418 421 1016 l6a 
478 474 . 77 6, 34 yp, 16b+ 10a 
664 669 1126 4 
729 T3177 Yugy, 13 6, 11 Vy @w(NO,)+4+ 10b 
797 794 99 yp 11 
891 886 107 yz 17b 
922 927 103 yp 17a 
980 976 119 yp 5 
para Bromonitrobenzene 
(b, species) 
102 91 23 Yr, 43 4, 38 yy 10b+ 16b+ 11 
276 251 38 yRy, 54 Yu 10b+ 11 
459 449 75 6, 33 yr, 16b+ 10b 
674 657 976, 16 yyy 4 
737 738 73 ym, 26 4, 10 yy w(NO,)+4 
951 951 102 yp 17b 
965 968 120 yz 5 
(a, species) 
405 440 1206 16a 
820 828 91 yp 10a 
843 839 100 yz 17a 


conditions with an average error of 8.7cm * between 
the observed and calculated frequencies. 


3 Results and Discussion 

The final set of force constants obtained, along with 
the initial set and dispersions, are presented in Table 1. 

The observed and calculated frequencies and the 
potential energy distributions (PED) in principal 
internal coordinates obtained as explained in an earlier 
paper’ are given in Table 2. Contributions below 10% 
are Table 2 contains vibrational 


= - not shown and 
7 ee inccts also. 
oe a += 


le 


pa } fi . 


al Assignments 
bles 1 and 2 it could be seen that: 
agreement between the observed and 


(ii) The final values of the force constants are close 
to their initial values as the corrections have never been 
very large, 

(111) The dispersions of the force constants are quite 
small, and 

(iv) The potential constants are comparable with 
their counterparts in related molecules. 

As such, all the requirements necessary for an 
acceptable force field are satisfied by the potential 
constants obtained for the out-of-plane fundamentals 
of mono bromonitrobenzenes. Hence, they are 
physically reasonable. Therefore, the assignment of 
various frequencies can be made on the basis of this 
normal coordinate analysis, rather than from 
qualitative considerations*”>. 

C —H out-of-plane bending vibrations—As identified 
from the phase relations and the PED presented in 
Table 2 (please see Ref. 7 for details), the bands at 985, 
775, 852 and 953cm~' are due to the C—H out-of- 
plane bending modes 5, 11, 17a and 17b respectively in 
the ortho compound. The corresponding vibrations in 
the meta compound are at 980, 797, 922 and 891cm =? 
whereas these are identified near 965, 820, 843 and 
951cm™' in the para isomer. Thus the mode 5 is the 
highest out-of-plane C—H bending vibration in the 
three molecules studied. This result agrees with that of 
Green and Harrison* for ortho and meta compounds 
but differs in the case of para isomer. The lowest out-of- 
plane bending mode is found to be that in which the 
corresponding internal coordinates are either 
increasing or decreasing simultaneously. It is 
established that frequency of mode 17b is greater than 
that of 17a in ortho and para compounds whereas the 
reverse is true for the meta compound. This is in perfect 
agreement with the expectations of Green and 
Harrison*. Further, the present computations reveal, 
all the C—H_ out-of-plane bending modes to be pure 
as is evident from Table 2. 

C—X out-of-plane bending vibrations—These are C 
—Br and C—N out-of-plane bending modes in the 
present set of molecules. It is seen from Table 2 that 
the bands at 132 and 199cm' in the ortho compound, 
the fundamentals near 147 and 193cm~'! in the meta 
compound and the absorptions around 102 and 
276cm~' in the para compound exhibit the nature 
expected of C—Br and C—N out-of-plane bending 
modes. These two modes mix with each other. It is 
interesting to note that the lower of the two modes has 
a large contribution of PED from the ring torsion 16b 
in the ortho and para compounds whereas the higher of 
the two gets such contribution in the meta isomer. 

Ring torsional vibrations—The three ring torsional 
modes 4, 16a and 16b are identified at 690, 413 and 
462cm~' respectively in the ortho compound (for 
details please see Ref. 7). The corresponding vibrations 
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in the meta isomer appear near 664, 418 and 478 cm Fi 


whereas in the para compound they are around 674, 
405 and 459cm~!. Thus the order of the ring torsions 
should be 4 > 16b > 16a in these three molecules. This is 
against the expectations of Green and Harrison* for 
ortho and meta compounds, viz. the frequency of mode 
16a to be greater than that for 16b. Further, the modes 
4 and 16a are pure in the three molecules whereas 16b 
has considerable contribution from C—Br out-of- 
plane bending mode. 

Out-of-plane vibrations of the nitro group—The out- 
_ of-plane wagging mode of the nitro group is identified 
at 730, 729 and 737cm ' in the ortho, meta and para 
compounds respectively. Thus, Green and Harrison’s* 
assignment for this mode in the ortho compound 

_ should be corrected from 775cm~! to 730cm_'. As 
mentioned earlier, the band at 775 cm ' is due to mode 
11 corresponding to C—H out-of-plane bending 

ee le ee J F 
vibrations. It is important to note that this mode mixes 


_ molecules under consideration. A discussion of the 
; Bi , ‘si onal mode of the nitro group is not given as it 
cot ulc | not be included i in the calculations for reasons 


Fitba highest C—H out-of-plane bending — 
eas that in which the corresponding internal | 
increase or decrease simultaneously is the 7 Mu ura 

: allt the CoH # out-of-plane tending = = 


additional PED contribution from ie ring to sion 


Ae 


16b. 
(iii) The order of the ring torsions is re 16b> 16 


Modes 4 and 16a are pure whereas 16b— 


considerable mixing with out-of-plane oe mode 7 


of:C — Br. ce 
(iv) Out-of-plane wagging vibration of the nit TO 


extent in the three molecules. 
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Dielectric constant (é’) and dielectric loss (e”) have been determined for pure 2, 4-, 2, 6- and 3, 5-lutidines at frequencies 
720 kHz, 6.118, 9.8 and 19.22 GHz for four different temperatures. Dielectric relaxation time (t) and the distribution parameter 
(a) have been evaluated by Cole-Cole arc plots obtained at temperatures 20, 30, 40 and 50°C. It is observed that a decreases with 
the increase of temperature for 2, 4- and 2, 6- lutidines. On the contrary it increases for 3, 5-lutidine with temperature. Various 
thermodynamic parameters have been evaluated by considering the relaxation as a chemical rate process. 


1 Introduction 

The measurement of relaxation time at different 
temperatures enables us to calculate various 
thermodynamic parameters which give very useful 
information regarding the molecular structure. These 
investigations give an idea about intramolecular and 
intermolecular forces. The distribution parameter (&) 
lies between 0 and | and it represents the distribution 
of the relaxation time. The larger « is, the larger is the 
extent of distribution of relaxation time. 

The various symbols and letters used in this paper 
signify as follows: 
T =temperature in absolute scale 
e =dielectric constant 
é” =dielectric loss 
Tt =macroscopic relaxation time 
T) =microscopic relaxation time 
=Planck’s constant 
= Boltzmann constant 
=gas constant 
= viscosity 
=distribution parameter 
=Avogadro number 
= frequency 
AF =molar free energy of activation 
AH=molar enthalpy of activation 
S =molar entropy of activation 


When ¢ and n appear as subscripts, they stand for 
relaxation and viscosity processes respectively. 


Sz2zRs by > 


2 Experimental Details ; 
The dielectric constant (e’) and dielectric loss (e”) of 
the pure lutidine compounds have been measured at 3 


microwave frequencies, viz. 6.118, 9.8 and 19.22 GHz 
seca csideures, viz. 20, 30, 40 and 50°C. 


method! was used for the determination of €’ 
” The static dielectric constant was measured by 
sole meter at 720 kHz at the above mentioned 


» 


temperatures. The apparatus works on the principle of 
heterodyne-beat method described elsewhere?’*. The 
waveguide transmission technique was used at 
microwave frequencies because this is most convenient 
and a fairly accurate method for measuring dielectric 
constant and loss in the liquid samples. 

High-purity lutidine compounds, viz. 2, 4-, 2, 6- and 
3, 5-lutidines were imported from Koch-Light 
(England), BDH Chemicals (England) and E. Merck 
(West Germany) respectively and were used as such. 


3 Results and Discussion 

The measured values of €’ and €” are represented by 
Cole-Cole* arc plots which are shown in Fig. 1[(A)- 
(C)] for 2, 4-, 2, 6- and 3, 5- lutidines respectively. Here 
€,, is taken as the square of refractive index. The t and 
the values of « evaluated from these plots are presented 
in Table 1 for all the 3 compounds. Besides this, the 
relaxation time for lutidine solutions using carbon 
tetrachloride as solvent (as a typical example) was also 
evaluated by Gopala Krishna’s method? at 9.8 GHz. 
The relaxation time in solution (t,,)) has also been 
compared with (t) in Table 1. 

The values of to calculated using Debye’s, Powles 
and Fatuzzo and Mason’s” relations are compared 
with the t values at 20°C (Table 2). The explicit 


relations for the to are given as follows® ~* 
os | Sale 
Debye relation Tt) = tot? T 7) 
. 29 + &, 
Powles relation To = . (2) 
3&5 
Fatuzzo and Mason relation tom tot te . 3) 


2& +8, 


Viscosity (n) values for all the 3 compounds have 
also been measured at different temperatures with the 
help of Ostwald viscometer. These 7 values and the 
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Fig. 1—Cole-Cole arc plots for: (A), 2, 4-; (B), 2, 6- and (©), 3, 5- 
lutidine 


interpolated values of densities are presented in 
Table 3. 

Various thermodynamic parameters for the 
dielectric relaxation and viscosity processes have been 
evaluated using Eyring’s equation’. In this case the 
dielectric relaxation is considered as a rate process 
defined by the equation: 


t=(h/kT) exp(AF,/RT) ... (4) 
where AF,=AH,— TAS, 


h AH, —AS 
T =| — — : 
hus Tt (Er) =s°( Fe )eo(=A5) ‘ Jh) 
the AF,, AH, and AS, being the same as defined in 
Section 1. 

Our observations support the theory of rate process 
for a dielectric. For all these compounds, the plots of 
_ a. 1/T and log(tT) versus 1/T gave straight 
Inés as demanded by the rate the ; 
Pe 2 Ory equations 

The AH, has been evaluated 
plot of log(tT) versus 1/ 
calculated from the relation 


AF, =2.303RT log (7) 


by using the slope of the 
Ty and AF, has been 


.. 6 
340 
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The values of AF,, AH, and AS, so evaluated are also 
presented in Table 1 foi the lutidines. The molar 
volume V(=M/d) calculated from molecular weight 
and density has been compared in Table 3 with the 
hypothetical volume (Vp) given by Debye’s equation - 


Vp=tRT/3n O 


The AH, has been calculated using the slope of the 
plot of logy versus 1/7 and the AF, from the relation: 


AF,, =2.303 RT log(yV/Nh) . (8) 


The AS, has been calculated from the following 
relation 


AF, =AH,— TAS, mC) 


The various molar activation energy parameters so 
evaluated are given in Table 1 for the 3 lutidine 
compounds. 

The factor V* which gives the dipole rotation 
volume?® i.e. the effective molar volume available for 
dipole rotation process, has been calculated from the 
knowledge of molecular weight (M) and density (d). It 
is given by the relation V*=(M/d)* where x 
=(AH,)((AH,)= ratio of enthalpies for the two 
processes and the values so calculated are presented in 
Table 3. For all the 3 compounds, V is found to be 
much larger than Vp, calculated from Debye’s 
equation, Vp being only 23-24.4% of V, 25.84-29.3% of 
V and 22-22.7% of V respectively as shown in Table 3. 
This is in agreement with the results obtained by other 
authors'*. The discrepancy is due to the fact that the 
Debye’s equation does not account for the deviation of 
molecular shape from the spherical symmetry. 
Another reason is that the macroscopic viscosity is 
different from the unknown microscopic viscosity, 


8-10 


8-00 


LOG TT 


3-00 310 3:20 3-390 3.49 a 
4% 103(K) : | 
Fig. 2—Plot of logy (dashed lines) and log tT (continuous line 
versus I/T forthe 3 lutiine compounds (A), 2, 4(B), 2. & ended 
3, S400“ ore he 
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Table !—Dielectric Relaxation and Molar Activation Energy Parameters of the 3 Lutidines Studied 
((A), 2, 4-lutidine, (B), 2, 6-lutidine, (C), 3, 5-lutidine] 


20°C 30°C 
(A) (B) (C) (A) (B) (C) 
&Cole —Cole 0.45 0.27 0.43 0.43 0.23 0.48 
(t x 10''), sec 2.90 3.35 ~ 3.40 2.49 2.95 2.80 
(t,.; x 10"), sec 5.12 7.50 5.30 4.87 5.95 4.55 
AF,, kcal/mol 3.02 3.10 3.11 3.05 3.15 3.12 
AF,, kcal/mol 3.23 3.25 3.33 3.27 3,275 3.36 
AS,, cal/mol/K —3.75 —4.76 =212 ~ 3,73 —4.77 —2.08 
AS,, cal/mol/K —2.36 — 1.64 —2.22 —2:41 — 1.67 —2.24 
40°C 50°C 
a 0.41 0.19 0.53 0.39 0.155 0.58 
(t x 10"), sec 2.23 2.60 2.37 1.96 2.32 2.05 
(To: X 10"), sec 4.49 5.55 3.93 = =e = 
AF,, kcal/mol 3.08 3.19 3.14 3.13 3.24 3.16 
AF,, kcal/mol 3.30 3.29 3.39 3.34 3.32 3.42 
: AS,, cal/mol/K —3.71 —4.74 —2.08 —3.75 —4.75 —2.07 
‘ AS,, cal/mol/K —2.43 — 1.66 +227 —2:48 —1.70 —2.29 
7 (A) (B) (©) 
AH,, kcal/mol 1.92 1.705 2.49 
AH,, kcal/mol 2.54 2.77 2.68 
“ x=AH,/AH, 0.76 | 0.62 0.93 Se 
a Up(static), D 1.67 caer | 1.50 1.70 ; <2 
spt microwaves), D- 1.65 1.45 1.67 
ee? 252. 1.84 < 2.67 ; = ae : 


| 1] measure of the internal friction, used Table 2—Macroscopic and Microscopic Relaxation Times 


I sion tion, . (7). The ratio of 7, to in the 3 Lutidines Studied _ 
, is th sar eat = ah | (All values at 20°C) 
F be observ sony reeks that for all the 3 Compound. ate i. 
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Table 3-—Comparison of Molar V 
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olumes in the 3 Lutidines Studied 


@, Dipole rotation 
=i 10!') Vp V_/V,% po : 
Temp. (7) Viscosity(”) Density(@) Ne al volume(V*) 
C cP g/ml “ 
2, 4-Lutidine 
90 26.56 23.10 36.85 
> 0.931 115.10 2. 
. a 0,922 116.23 2.49 26.75 23.01 a 2 
0 0.676 0.913 117.40 2.23 28.63 24.40 
= 0.607 0.903 118.63 1.96 28.91 24.40 37.70 
2, 6-Lutidine 
30.00 25.84 19.10 
0.907 0.923 116.1 3.35 
e 0.780 0.910 117.8 2.95 31.80 27.00 19.24 
40 0.661 0.900 119.1 2.80 34.13 28.66 19.37 
50 0.583 0.880 121.8 2.32 35.63 29.25 19.64 
3, 5-Lutidine 
20 1.073 0.943 113.64 3.40 25.74 22.70 81.60 
30 0.912 0.932 115.00 2.80 25.80 22.43 82.50 
40 0.795 0.921 116.40 2.37 25.87 22.23 83.43 
50 0.700 0.910 117.80 2.05 26.22 22.26 84.40 


more uniform at higher temperatures. The increase in 
the thermal agitation favours predominance of those 
configurations which have more nearly equal free 
energies of activation. It may be said that in 2, 4- and 2, 
6-lutidines interaction between methyl groups and 
nitrogen atom is quite significant. Therefore various 
segments of the molecule relax together. Thus 
cooperative effect of relaxation occurs in both cases, 
and the effect increases with temperature. These results 
are in agreement with those reported in 
literature’? ~!9, 

Itis further observed from Table | that in the case of 
3, 5-lutidine & is increasing with temperature, contrary 
to the other two lutidines. In the case of 3, 5-lutidine 
the cooperative effect of relaxation seems to decrease 
with increase of temperature. Similar observations 
have also been reported in literature by other 
workers'*''®. This result can be explained in the 
following manner: 

In the case of 3, 5-lutidine two methyl groups are 
attached in 3, 5 positions with respect to the nitrogen 
atom, i.e. they are comparatively away from nitrogen 
atom and there is a weak interaction between them. It 
Is expected that when this compound is excited at 
higher temperature, thermal fluctuations (vibrations) 
are more dominant as compared to this weak 
interaction. Hence two methyl groups and other 
segments of the molecules are influenced by the 
thermal vibrations and they relax independently at 
higher temperature. In other words, as the temperature 
increases, various segments of the molecule have more 
thermal energy and they are now capable of relaxing 
independently. 

The molecule of 3, 5-lutidine may be visualized as 
relaxing due to: (1) the rotation of one methyl group (2) 


the rotation of another methyl group (3) the rotation of 
H of CH and (4) the overall molecular rotation. Thus at 
higher temperatures this molecule has four very 
distinct ways of relaxation which accounts for the 
increase of « with temperature. 

It is evident from Table 3 that V* available for dipole 
rotation process increases with temperature in each 
case viz. 36.85 to 37.70 cc, 19.10 to 19.64-cc and 81.60 
to 84.40 cc for2, 4-, 2, 6- and 3, 5-lutidines when the 
temperature is increased from 20 to 50°C. Evidently, 
the effective molar volume available for dipole 
rotation process is largest for 3, 5-lutidine and smallest 
for 2, 6-lutidine. Consequently, the freedom of dipole 
rotation is greatest in 3, 5-lutidine which gives rise to 
the increase of « with temperature. 

The process of dielectric relaxation and the process 
of viscous flow are analogous because in both the cases 
the molecule requires an activation energy to 
surmount the hindering potential energy barrier. A 
comparison of the two energies using Table 1 reveals 
that for the compounds studied the AF,, is greater 
than AF,. This is expected because the viscous 
flow involves the translational and rotational motions 
of the molecule whereas the relaxation process 
involves only the rotational motion. Both AF, and AF, 
show a gradual increase with temperature. It is 
observed that for all the 3 compounds the molar 
enthalpies of activation AH are smaller than the 
corresponding free energies of activation and hence the 
most probable entropies of activation AS are found to 
be negative. All these results are in agreement with 
those reported by other workers!? 15. 4 to 
Branin and Smyth'* a negative entropy of activation 
indicates that there are fewer configurations possible 
in the activated state and for these configurations the 

¢ “ 
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activated state is more ordered than the normal state. 
This may be interpreted as indicating the presence of 
the cooperative orientations of the molecules resulting 
from the steric forces or from the dipole-dipole 
interactions with the dipoles being more nearly aligned 
with each other in activated state. 

It is also found that for 2, 4- and 2, 6-lutidines the 
change in entropy due to relaxation process is greater 
than that due to viscosity process. Further, this change 
in entropy of activation is greater for 2, 6-lutidine than 
that for 2, 4-lutidine. This suggests that in 2, 6-lutidine 
an inductive effect may be more significant because of 
close proximity of two methyl groups to the nitrogen 
atom. This inductive effect causes more intramolecular 
polarization, mainly orientational polarization and 
hence there is greater change in entropy, the latter 
being a measure of maximum disorder. For 3, 5- 
lutidine the change in entropy of activation due to 
viscosity process is found to be greater than that for 
relaxation process as AS, is more negative than AS,. 

The ratio of enthalpies of activation for the two 
processes, viz. (AH,/AH,) is found to be 0.76, 0.62 and 
0.93 for 2, 4-, 2, 6- and 3, 5-lutidines (Table 1). 
Therefore, according to the classification rules 
suggested by Krishnaji and Man Singh’”, it is expected 
that all of them will belong to class I liquids in which 
the freedom of dipole rotation is immediately lost on 
solidification because dipole freezing will occur on 
solidification. Consequently, there will be a sudden 
drop of dielectric constant for these compounds at 
their freezing point. 

The values of dipole moments (ps) have been 
calculated by Debye’s equation at static and 
microwave frequencies and by Onsager’s formula at 
static frequencies as presented in Table |. These 
formulae are as follows: 


(Debye) p45 (Static) 
9kTM E —1 


———_—_- 


te hae 
Eg +2 ; 


= 


~ 4nNd 


(Onsager) 46 (Static) 


_9kTM (€9 —€~)(2&9 + Ex) .. CD) 
~ 4nNd El Eq + 2) 


(Debye) “5 (Microwaves) 


_27kTM ¢' 1+m?*r2 (12 
~ 4nNd (é9+2)(e,+2) wrt =< 


where €,, is taken as the square of the refractive index. 

As expected, the values of fips calculated by Debye’s 
equation are found to be lower than those calculated 
from Onsager’s formula. This is because in the 
Onsager’s formula due consideration has been given to 
the presence of the internal field. The values calculated 
from Onsager’s expression are in good agreement with 
those obtained for solutions. 
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Dielectric measurements have been made on some rigid polar molecules in dilute ee of ee ee 
mixture at 9.8 GHz microwave frequency and over a temperature range of cig K. : fe tae A le 
resolved using the dielectric data by applying the methods of Higasi, and Higast, Koga an is Fi a fake 2 The 
the distribution parameter obtained for these systems indicates the possibility of hindered motion cal ve ; . Z 
flexible behaviour of these solutes has been interpreted in terms of recoiling and consequent trapping of t e Ye) ute mo ecu 3 
between the solvent holes. The hindered motion observed in highly viscous medium in the case of the investigated rigi 
molecules exhibits good agreement with the earlier study of G A Balogun and C W N Cumper [J Chem Soc Faraday Trans I 


(GB), 69 (1973) 1172]. 


1 Introduction 

The effect of macroscopic viscosity on the relaxation 
time of rigid polar molecules has been investigated by 
several workers' °. Kalman and Smyth* observed 
that the viscous solvents tend to increase the 
separation between the two relaxation processes and 
the distribution parameter is considerably increased in 
highly viscous medium. Balogun and Cumper® 
Suggested that the dielectric relaxation in polar 
molecules occurs by their rotational diffusion and its 
rate depends upon the magnitude of intermolecular 
forces present in the viscous solution. They studied 
methyl iodide, cyanocyclopropane, chlorobenzene, 1- 
chloronaphthalene, etc. in different mixtures of 
paraffin oil and benzene of varying viscosity and 
observed that even rigid solute molecules in a highly 
viscous medium exhibit a small unsymmetrical 
distribution of relaxation times and that a Debye 
treatment based on single relaxation time would no 
longer be strictly adequate. Keeping this in view, the 
dielectric measurements on four rigid polar molecules 
(pyridine, quinoline, bromobenzene and iodobenzene) 
in six solvents of varying viscosities were made by us at 
@ microwave frequency of 9.8 GHz and at a static 
frequency in the temperature range 298-322 K. The six 
solvents were Prepared by adding increasing amounts 
of paraffin oil to benzene. 

The distribution parameter « and the most probable 
relaxation time Ton were evaluated using the 
hour mio seat Sic 

;, appreciably high values 
of “* in highly viscous Media it was considered 
ee - resolve the relaxation mechanism into 

(1) second representing the molecular 
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rotation T,,) using Higasi er al. method’. The average 
relaxation time to, has been calculated using the 
equation T,o)=[T 1) T,2)]"/*. The details of evaluating 
these parameters have already been discussed in our 
earlier papers*’’. The observed results are discussed in 
the light of Balogun and Cumper’s model® of hole 
creation and the subsequent trapping of the rigid 
solute molecules. 


2 Experimental Details 

The following six solvents of varying viscosities were 
taken and the solutions were prepared by changing the 
amounts of the solute dissolved: (1) 100% benzene; 
(2) 80% benzene + 20% paraffin oil: (3) 60% benzene 
+ 40% paraffin oil; (4) 40% benzene + 60% paraffin oil; 
(5) 20% benzene + 80% paraffin oil; (6) 100% paraffin 
oil. The experimental details of the measurements and 
evaluation of dielectric parameters have been given in 
our earlier papers® ~!°, 

The chemicals pyridine, quinoline, bromobenzene 
and iodobenzene were obtained from M/s British Drug 
House, England. The medicinal paraffin oil was 
obtained from Liberty Pharmaceutical Ltd, Bombay; 
the other solvent benzene (AR grade), obtained from 


M/s British Drug House, England, was distilled twice 
before use. 


3 Results and Discussion 

The relaxation times, evaluated using the Higasi® 
and Higasi er al. method’, and the distribution 
parameters of all the four compounds are listed in 
Table | along with the macroscopic viscosities of the 
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Distribution Parameter («) and Macroscopic Viscosity (Ys in mP) 


Table 1—Relaxation Times (rt in ps), 


Temperature TH T2) To) Ton) x Ns 
K 
q I. PYRIDINE 
| 100% Benzene 
314 5.4 7.9 6.3 5.1 0.08 4.8 
322 4.3 $7 4.9 4.1 0.06 3.3 
80°, Benzene + 20% paraffin oil 
298 8.2 11.7 9.8 8.4 0.11 — 
; 306 7.2 8.4 78 e213 0.05 10.8 
P 314 5.3 71 6.1 5.5 0.07 8.8 
322 5.2 5.8 “bigs 5.3 0.01 6.6 
60%, Benzene + 40% paraffin oil 
298 9.9 20.6 14.3 11.3 0.23 — 
306 9.6 15.6 12.2 9.0 0.15 16.3 
314 8.1 10.2 9.1 8.1 0.07 12.0 
322 a 8.3 8.0 5.4 0.03 8.9 
: 40% Benzene + 60", paraffin oil 
298 6.7 31.1 Mae aS 0.45 — 
306 6.5 28 pee See “a2 0.42 28.7 
314 = (25.7 11.5 83 0.40 20.7 


poe ae eee: 


20”, Benzene +80", paraffin oil 
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and Macroscopic Viscosity (7s in mP) | 


Table |—Relaxation Times (t in ps), Distribution Parameter (x) 
Temperature tp Te To) T.OH) x Ns 


20°,, Benzene + 80", paraffin oil 


34.4 
34.2 
32.5 
21.4 


41.2 
39.9 
36.5 
31.4 


III. BROMOBENZENE 


11.6 
10.7 
9.8 


17.8 
17.3 
15.5 
12.0 


100°, Paraffin oil 


18.3 
17 
14.8 
13.0 


100°,, Benzene 
10.3 


9.2 
8.9 


22.8 
19.4 
18.7 

8.0 


23.7 
21.3 
16.0 
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9.6 
8.0 
6.8 


80°, Benzene +20° » Paraffin oil © 


16.7 
15.9 
14.6 
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SINGH et al.: DIELECTRIC RELAXATION & MOLECULAR MOTION IN HIGHLY VISCOUS MEDIUM 


Table |— 


Relaxation Times (t in ps), Distribution Parameter (#) and Macroscopic Viscosity (7s in mP)—Contd 
Temperature wtp T2) To) OH) m Ns 
K 
60% Benzene + 40% paraffin oil 
298 10.7 28.9 17.6 19.2 0.31 
306 10.6 27.8 FT 19.1 0.29 16.3 
314 9.6 22.8 14.8 13.5 0.27 12.0 
322 8.8 16.7 12.1 12.3 0.15 8.9 
40%, Benzene + 60% paraffin oil 
298 11.3 32.8 19.2 21.5 0.32 — 
306 11.2 29.2 18.1 20.4 0.30 28.7 
314 11.1 26.5 17.2 18.0 0.28 20.7 
322 13.8 21.8 17.4 17.9 0.15 15.9 
20%, Benzene + 80% paraffin oil 
298 11.7 38.6 21.3 21.7 0.31 = 
306 9.4 30.1 16.8 21:5 0.43 104.6 
314 9.2 29.0 16.3 16.5 0.35 95.1 
322 9.4 26.0 15.6 14.9 0.32 78.1 
100% Paraffin oil 
298 8.2 41.2 18.4 24.2 0.47 — 
306 7.9 41.1 18.2 22.8 0.48 483.4 
314 7.6 38.4 17.3 19.1 0.47 384.6 
322 7.6 36.9 16.7 17.9 0.46 240.3 


viscosity. As may be seen from Table | for the pyridine 
molecule, t,o} (at 314 K) varies from 5.1 to 10.8 ps in the 
range of medium viscosity of 4.78 to 384.6 mP whereas 
To) increases from 6.3 to 14.1 ps in the same viscosity 
range. It is also observed that the two relaxation 
parameters t,o, and tio) are comparable. A similar 
behaviour is observed in the case of the other three 
molecules investigated here. 

Since the dielectric data yield substantial values of 
distribution parameter in the range 0.01-0.47, with 
increasing viscosity, the process was further resolved 
into two separate relaxation times 7,,) and 7,2) using 
Higasi et al. method’. The observed high value of « 
indicates the possible flexible behaviour of these rigid 
molecular systems. 

A comparison of T,,, and 7,,) for the pyridine system 
exhibits that the two parameters are nearly equal in the 
dilute solutions of benzene. With increasing medium 
viscosity, 1,2) increases systematically giving a 
maximum value for the solutions having 100% paraffin 
gil. A similar behaviour is exhibited by the other three 
rigid molecules under investigation. On the other 
hand, the variations observed in the case of T1) values 
are small and are accompanied by an irregular 
behaviour with increasing viscosity of the medium. As 
suggested by Balogun and Cumper, the solvent 
molecules in a highly viscous system form a type of 
cage around solute molecules and for the dipolar 
rotation to occur in the solution, the solute molecules 
have to pass through various hindrances created due to 
ge. The varying degree of freedom 


experienced by the solute molecules would result in 
different values of the relaxation times. In the present 
case of rigid molecules, t,,, and 7,2) are such that T,;) 
would possibly represent the relaxation time of the 
part of the segment of the solute molecule recoiling 
around the solvent molecule and 1,,), the relaxation 
time of the solute molecule as a whole performing the 
dipolar rotation. The irregular behaviour of 7,,), with 
variation of viscosities, may be ascribed to the small 
hindrances produced from the segmental motion in the 
solvent cage. On the other hand, a systematic increase 
of the t,) values with increase in viscosity would be 
expected due to the dipolar rotation of the entire 
molecule as a unit. 

Another rigid system, viz. quinoline, exhibits a 
similar behaviour. It may be seen from Table | that the 
T;,) and t,2, values for the temperatures investigated are 
not substantially different in benzene solution. Again, 
Tj) gives an irregular variation with increasing 
viscosity whereas T,) increases systematically giving 
values of 16.0-41.2 ps in the viscosity range 5.6-240.3 
mP. The other two systems, viz. bromobenzene and 
iodobenzene, also exhibit a similar behaviour. 

These investigations indicate that the distribution in 
the relaxation times occurs even in the case of rigid 
systems in highly viscous media. Here the systems 
taken were both spherical and slightly elongated types 
and the phenomenon seems to hold true equally in 
both the systems. As has been mentioned earlier, the 
behaviour is perhaps governed by the concept of cage 
formation of the solvent systems and consequent 
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trapping of the:solute molecules. The strength of the 
hindrance induced depends upon the number of such 
cages created which varies with solvent viscosity and 
which would effectively contribute towards the barrier 
provided for both the molecular motion and the 


segmental flexibility. 


4 Conclusion 
All these four rigid polar molecules in the high 
viscous medium showed sufficient amount of hindered 
motion. This hindrance is perhaps due to diffusion, 
recoiling and trapping of solute molecules in the 
_ solvent molecules. The observed high values of the 
_ distribution parameter («) would arise from these 
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B V GURUNADHA RAO, V P BHATNAGAR* & S B AGARWALt 
Department of Physics. Delhi College of Engineering, Delhi 110006 


Received 2 January 1982; revised received 21 May 1983 


Mutarotation of fresh and equilibrium solutions of glucose have been studied under irradiation of ultrasonic waves of 
1.5 MHz. The process of mutarotation has been found to be accelerated under ultrasonic exposure. Further, a marked decrease 
in the optical rotation was noted in case of equilibrium solution. This effect increased with the exposure time to ultrasonic waves 
and attained a saturation. The time-constant by which the glucose solution attained its new equilibrium value was found to 
depend on ultrasonic power. The observed effect indicates that the structural state of the glucose can be changed by ultrasonic 


irradiation. 


1 Introduction 

It is a well known fact that when glucose is dissolved 
in water, the optical rotatory power of the solution 
decreases on standing until finally it reaches a constant 
value, a-phenomenon which came to be known as 
mutarotation. This change in optical rotation has been 
attributed to the change in molecular structure of the 
glucose solution. This explanation received a striking 
confirmation in the discovery of two forms of this sugar 
due to asymmetry in the reducing carbon atom 


-mumber-! of the cyclic sugar. The alpha and beta 


modifications of the various ring forms and the open 


<. | exist insolution in a dynami equilibrium!. The optical 
rotations of these substances are not alike and _ 


iently when a sugar is dissolved in water, the 
a changes as equilibrium is established 
‘ee ae The dug tag? of 


sugar and may facilitate the separation of different 
modifications of the sugar. d-Glucose crystallises from 
pyridine solution in the form of a pyridine compound 
which contains the beta pyranose modification. But it 
crystallises from water solution at ordinary tem- 
peratures in the form of a hydrate which contains 
the alpha pyranose modification. The mutarotational 
equilibrium thus changes with chemical process. 
Sometime back, two of the present authors have 
reported that the optical rotation of glucose is changed 
under physical processes 
ultrasonics?*. The experimental data in earlier 
experiments were obtained by visual polarimetric 
techniques, without due considerations being given for 


few aun it is essential to use a solution tha 
attained equilibrium so as to get meanir 

Also it was of interest to study the time v. 
ho cadet Pome eeerte 


“nga mi sive ad ae Pate 
= ae (eet A using: an» ‘2 


sar polatir onic irradiation = ae Dees 


ar ar ‘ay are pre eho int 


ae 


like propagation of 


the age eS solution. ‘Since pe gianna of glucose 
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The variation in optical rotation of the 36:64 glucose 
solution with exposure time was studied for different 
ultrasonic intensities. A family of curves so obtained is 
presented in Fig. 3. ere 

The curves show that-a new equilibrium is reached in 
an exponential manner according to the following 


equation: 


Fig. |—Experimental set-up to study mutarotation under ultrasonic acstha Rees 
irradiation [1 Mercury source, 2Green light filter, 3 Polariser, 
4 Analyser, SExperimental cell, 6 Thermometer, 7Transducer, where t is a time-constant and the normal value 2% 


8 Photomultiplier, 9 Microvoltmeter, 10CRO and 11 Ultrasonic correspon ding to.t= eh given by 
power generator ] 


Ps % =(A+B) 
transversely the acoustically disturbed medium was a5 . 
observed using nicols fitted in circular scales. A and the new equilibrium value «,, corresponding to 
photomultiplier and a microvoltmeter were used to [t=© is given by 
read the microscopic changes in the optical rotation. saad 
For this purpose, a graph between microvoltmeter = Sas a 
reading (x) and rotation angle (0) of the analyser was From the curves in Fig.3 it is clear that the new 
plotted and a point giving maximum sensitivity, ie. equilibrium value attained under different intensities is 
(dx/d6),,,. was selected as the operating point to study the same for a particular specimen, however, the half- 


the minute changes in the optical rotation. A sensitive N15 
thermometer was kept in the experimental cell to read H,0 
the temperature variations of the solution during the wo 


ultrasonic irradiation. A oscilloscope was also used to 


check the frequency and to measure the rf voltage o 
across the crystal. ae 
oO 
a 225 
3 Results and Discussion K 60 Pa 
The experiment was done in two parts. In the first : —e 
part, the fresh solution of glucose was exposed to 4% 
ultrasonics and in the second part, the experiment was 40 
done on that solution of glucose which after natural 
process of mutarotation attained the equilibrium. We 
shall call this solution as 36:64 glucose as it is well 0 ry 7s nt 5 2 re 
known that it contains 36 % of -glucose and 64 * of B- TIME , min : 


' Fig. 2—Optical rotation of fresh glucose as a function of exposure 
The first reading of optical rotation could not be time and ultrasonic power as €xpressed in terms of crystal rf voltage 


taken earlier than 3 min with fresh solution. On 
exposing the solution to ultrasonic waves, the value of 
optical rotation decreased with a faster rate than 
otherwise. In this respect the action of ultrasonic waves 
was similar to that of the catalytic action of adding 
dilute acids and alkalies. The effect of ultrasonic waves 
was studied at various intensities and it was found that 
the acceleration was greater at higher intensities 
(Fig. 2). 

With equilibrium solution, a marked decrease in the 
value of optical rotation was observed on ultrasonic 


and attained a saturation value, which was m hl 
than the value attained in the natural Process of 


the earlier communication no consideration given 
to the ageing of solution, with the adie an 


| 
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Fi ig. 4— Variation of time-constant with ultrasonic power as 


expressed in terms of crystal rf voltage 


od time is different for different intensities. The 
period times were determined from the slope of the 
line obtained on plotting log (x— A) agamst ¢. 
srimental study also reveals that lower the 


c Pesca, hee is the ee time, i.e. 


a “rotation is” an iitioation of the 
solution. Bd a apnoea 


> MUTAROTATION OF GLUCOSE UNDER ULTRASONIC IRRADIATION 


than a decrease in the optical rotation as has been 
observed. 

Petrenko® theoretically explained the ultrasonically 
induced optical rotation as a physical phenomenon, in 
certain substances. We feel that this explanation is not 
valid in the case of glucose solution because the effect 
persists even after the withdrawal of ultrasonic 
irradiation. 

Since the action of ultrasonic waves on freshly 
prepared solution is similar to the catalytic action of 
acids and bases which accelerate the process through 
the addition of H* and OH ‘ ions, it appears that sono- 
ionization of water of glucose solution is responsible 
for the ultrasonic action of speeding up the process of 
conversion of «-glucose into f-glucose to attain the 
equilibrium condition quicker. Also when equilibrium 
solution is irradiated by ultrasonic waves, the process 
of conversion of a-glucose into f-glucose continues 
further until the whole of the a-form of glucose 
transforms into the f-form. 
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Excess volumes (V*) and deviations in isentropic compressibilities 
(K,) have been determined at 30°C for the binary systems of 
acetonitrile with n-pentanol, n-hexanol, n-heptanol and n-octanol. 
Excess volumes have been measured directly by the dilatometric 
method and isentropic compressibilities have been computed from 
experimentally determined densities and velocities of sound. The 
positive values of V* and K, in all the systems are attributed to weak 
complex formation and predominating dissociation of associated 
alcohols by acetonitrile. Values of V’ and K, increase with an 
increase in chain length of alcohol. 


The nature and degree of interactions that are present 
in the binary mixtures of acetonitrile with n-propanol, 
_tpropanol, n-butanol, i-butanol and cyclohexanol 
have been Pee alitatively in terms of excess 
vo aa eae nae The Present 


d n-octanol. These alco 2 have been 


~<a eaanalabensatag but for the — 
- as ‘oyapahies with n-pentanol, n-hexanol, n- 


0-2 0-6 > 0-8 7-0 


hexanol; 3, pales and 4m — sae ; 
Allee 


The isentropic comprestibities () | 
calculated from the equation — ee 
y a aT eae $y aa 


= pares | 
k, ue “p ae rte bas Be 


whee. u 1 and f denote sound 
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Table 1—Values of Parameters bo, b; and b, in Eq.(4) 
System b, b, 
Acetonitrile + n-pentanol 109 -—39 —7I 
Acetonitrile + n-hexanol 130 =-5 —44 
Acetonitrile + n-heptanol 136 18 —34 
Acetonitrile + n-octanol ISS —-ii 1 


Standard deviation o(K,)=1 in each system. 


chain length. The strength of hydrogen bond in the 
aggregates and also the complex formation increases 
with decrease in chain length. Consequently, when 
homologous series of alcohols are mixed with a 
common solvent like acetonitrile the expansion in 
volume would be minimum in the case of lower alcohol 
and maximum in respect of higher alcohol. Similarly, 
the contraction would be maximum in the case of lower 
alcohol. Hence the observed values are in order. The 
order for values of K, may also be attributed to Ae 
same factors. | 

One of the authors (GDR) is ‘grateful to the 
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Amongst different approximate methods for the determination of 
interface state density at the Si-SiO, interface, the Hill and Coleman 
method is found to be the most suitable one. Results obtained from 
various samples having surface state density of varying magnitudes 
computed by three different methods are compared to support this 
conclusion. 


Amongst different methods proposed in the literature 
for the analysis of interface state density! * at the Si- 
SiO, interface, MOS parallel conductance technique 
developed by Nicollian and Goetzberger® is the most 
accurate one. Since the analysis of the parallel 
conductance data is quite elaborate and difficult, it is 
unsuitable for quick estimation of surface state density 
of a system. Several methods for the quick estimation 
of surface state density at the interface have been 
Suggested by previous workers®. These methods are 
approximate and give results which vary in magnitude 
by large amounts, say by factors of 2 to 4. Recently Hill 
and Coleman’ have proposed a single high-frequency 
approximation technique which is not only simple and 
quick in analysis but also provides results within 30 % of 
the most accurate one. In this note, results obtained by 
three different methods are compared. Various samples 
having surface state density of varying magnitudes 


Table 1— 
Sample Doping Crystal Oxide 
No. density orienta- thickness 
(cm~ +) tion (A) 

MOS, aS X 19? «111» 1000 
MOS, 1.8 x 10'® «111> 500 
MOS, 2.0 x 10'® €111> 1100 
MOS, 1.5 = 10'5 €1005 800 


**For PSG (phospho-silicate glass) de 
passed for 5 min. This was followed by 
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Position, a mixture of N, (1.5 x 10 
N, treatment at 1000°C for 30 min. 


were analyzed. However, representative results of only 
four samples are given in this work. 

Experiments were carried out on Si wafers of n-type 
having different resistivities and crystal orientations. 
They were chemically polished on one side and treated 
with hot sulphuric acid and dilute hydrofluoric acid (10 
parts water and | part HF). Finally they were cleaned 
with high resistivity deionized water. Diffusion of n* 
was made on the unpolished surface. The polished 
surface was protected by thick oxide previously grown 
on it. After cleaning in vapour degreaser the wafers 
were loaded in oxidation furnace for growing oxide on 
the polished surface under the appropriate ambient 
and temperature conditions. Aluminium dots were 
deposited on the oxide using vacuum vapour- 
depositing technique. Oxide thickness was determined 
by ellipsometer and checked from the magnitude of the 
capacitance of the structure under strong accumu- 
lation. A summary of the oxidation and annealing 
conditions is given in Table 1. 

Capacitance and equivalent parallel conductance of 
the MOS capacitors were measured with the help of 
H.P. multifrequency LCR meter model 4274A at 
several frequencies from 100Hz to 100kHz. The ac 
signal amplitude was not allowed to exceed 25mV so 
that the non-linear effects were negligible. Fast surface 
state density (N,,) at the Si/SiO, interface was obtained 
using Hill and Coleman’ single-frequency approxi- 
mation formula 


N = oe Gr, enaxs 0 
= q-A (Gin, max/@C x) +. =Cy/Co4 


where Gi max iS the measured conductance maximum 


Particulars of Si Oxidation and Annealing Conditions for Various Samples 


Field plate Annealing and oxidation 


area, A ___ processes used 
(cm?) 
5.1.x 10°* Dry oxidation at 1000°C 
491x10°° Dry oxidation at 1000°C 
+annealing for 20 min in N, 
ambient at same temperature 
3.02 x 10°? 0,4+5% gaseous HCl at 1000°C 
followed by annealing in N, 
ambient for 30 min at 1000°C and 
PSG** deposition 
6.36 x 1073 


O,+5% gaseous HCl at 1150°C 
followed by annealing in N, 
ambient for 40 min at 1000°C 
+PSG deposition 


ce per min) +O, (30 ce per min) + POCI, (30 cc per min) was 


NOTES 


in G-V plot, C,, the corresponding capacitance value 
on the C-V plot, A the metal plate area, w the angular 
frequency used, q the electronic charge and C,, the 
capacitance of the structure under strong 
accumulation. 

Traditional Nicollian and Goetzberger> analysis 
corresponding to any desired value of surface potential 
‘¥, was used for determining the surface state density at 
the Si-SiO, interface. For this, one usually determines 
C-V and G-V plots at several modulation frequencies. 
From these C-V and G-V plots and for any desired ,- 
value, G,/@ are calculated for the various signal 
frequencies used*. Next, a plot of G,/w versus @ is 
made. From the peak value of G,/w versus @ curve, one 
finds the interface state density as 

Ci 


~ cans A 


where C,, is the surface state capacitance, given by 


c(2). 


The surface state density was also determined using 
Castagne and Vapaille method’. Specesaine: to this 
eativod, N,, is given as 


, i Ean fu (= Cir/Cox Cue/Cox ) 
GANA —Cre/Cx)  (—Cur/Cos) 


Cir is the low-frequency capacitance and Cur 
1e h quency capacitance. — 
Fou r different samples having different diate state 


nt methods of Castagne-Vapaille, 
-Goetzberge and Hill-Coleman. The results 
apared | in nates Det ofthe eprinental 


at the Si-SiO, interface were analyzed by the | 


conductance curves for all the samples pertain to the 
depletion region. In Fig.2, G,/wC,, versus f for the 
structure at the gate bias —2.4 V is plotted. The bias 
corresponds to the conductance peak in the G,, versus 
Vg curve (Fig.l) of the device at 100kHz signal 
frequency. 

The single frequency approximation method 
developed by Hill and Coleman’ is found to be more 
suitable than any other approximation method. It may 
be pointed out that Hill and Coleman approximation 
method uses single time-constant for the surface states 
whereas in Nicollian and Goetzberger method 
provision for spatial fluctuations of surface charge can 
also be made, leading to the inclusion of time-constant 
dispersion effects in their framework. Nevertheless, 
Hill and Coleman’s method, being simpler and 
providing reasonably accurate results, is quite suitable 
especially for systems having low surface state density. 
Further, the surface potential dependence of time 
constant should make only a small contribution to the 
time constant, especially if one is confined to a small 


- 
s 
°o 


CAPACITANCE (Cm), pF 
ay 


a 
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High Field Single Injection Current in transition plane x, is shifted towards cathode when the 
Solid State Diode with Shallow Traps electric field increases. The low field region (0 < x < 
by Exact Method x.) lies adjacent to cathode and the high field region (x, 


S x S L)exists adjacent to anode. At sufficiently high 
A S VERMA, Y K SHARMA*, R S AGRAWAL & CS _ fields, the low field region is negligibly small and the 


AGRAWAL insulator is working under high field conditions. 
Department of Physics, R B S College, Agra 282003 The exact current-voltage characteristic may be 
Received 10 May 1982; revised received 22 September 1982 conveniently obtained with the help of the following 


dimensionless variables: 
An exact method is used to evaluate the current-voltage 1/2 
characteristic of a single injection solid state diode with shallow traps — No a NoMolLE.A(x)] 
Operating in a high field regime. A set of curves are drawn for various ohn J ; 
values of the parameter characterizing the shallow traps. 2 eotae 
Nou EL?x e NoHoE. Vix) 


Injection of current in an insulator has been widely ~ os ed O Oed? E(x) aEe 
studied in the presence of various physical parameters tee a 

of interest operating under high field conditions! ~°. The critical dimensionless variables at the transition - 
The mobility is field dependent at high fields. The low Plane x. are given from Eq. (5) 

field region is always present inside the insulator near eNoltoE- e? No UoX- 


cathode? which is generally neglected by the workersin  “«~ pag at er epars # 
the field. Here, the low field region is also included to oe 
study the current-voltage characteristic of the v= e Nolo ene ef : if ) 
insulator with shallow traps at high fields. OeJ | eee 
<< agpegerate insulator containing a significant — Consider the injection level at which ase on Hell Zs 
z of shallow traps which generate the and high field regions are present inside the insulator, 
Btration 2) nisouetg generated free carriers No. From Eqs(1)-6), the dimensionless os at pocsye 
lations characterizing the current flow derived as | ee 
Taw. high fields are given. si : ee a 


a Saigo ntee 
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Fig. |—Exact current-voltage characteristic of an insulator with 
shallow traps at high fields for 92=10~°, 10 ~* and 107 


scale in Fig. 1 with the help of Eqs. (7){9) where the 
critical dimensionless variables u,, w, and v, are 
neglected by assuming that the injection level of 
current is sufficiently high. The dimensionless values 
1{u,w,)''? and v,/w2 are obtained by the numerical 
tabulation of w, and v, against u,. The range of 
variable u, lies between 0 and 1. Other physical 
parameters are selected for cadmium sulphide at room 
temperature: 


Ng=10'°m ~*, L=10->m, e=10-!°F/m 
Mo =2 x 10 -7m2/V-s, E,=10’ V/m 


and 6=10~°, 10 * and 10 2 vee (LQ) 


The current-voltage characteristics of three samples 
with different values of @ are drawn in Fig. 1. It is 
shown that all the curves initially obey Ohm’s law (J 
x V) and they finally merge into a non-linear current- 
voltage characteristic. The effect of traps appears on 
the non-linear portion of the continuous curve, i.e. ata 
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higher injection level. The location of current-voltage 
characteristic is dependent on 9. The characteristic 
located on the right side of Fig. | is for lower values of 
s In the present problem, the presence of shallow traps 
is taken into consideration. The relation between the 
trapped electron concentration (n,) and characteristic 
constant @ is given by* 


(n/n,) = 9 wan (11) 


where n is the free electron concentration. In the 
absence of shallow traps, 0 is unity in Eqs (2) and (9). In 
general, the functional dependence of n, on n for any 
trap configuration is given by the simple form*: 


it a 


n,ocn? 


which is valid for the limits 0<p <1. Let us consider a 
shift in the quasi Fermi level F > F + oF. There will be 
a maximum change in trap occupancy if traps are 
shallow, n,  n, so that p= 1; and the traps are deep for 
minimum change, i.e. no change at all, in which case 7, 
=constant, so that p=0. 

If there are several sets of shallow trap levels in a 
given sample, the set with the smallest 0 dominates the 
current flow and this value of 6 satisfies Eq. (9). It may 
be noted that the traps are shallow so long as the quasi 
Fermi level lies below the energy level of traps. 
However, the quasi Fermi level keeps rising with 
injection level, i.e. with applied voltage. The effect of 
shallow traps on current-voltage characteristic is 
observed through the parameter 0. 
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The transmission of sonics through an ionized medium has been 
numerically evaluated with the help of a computer. The analysis 
shows that if a source of noise is surrounded by a plasma barrier it 
can effectively shield most of the noise spectrum. 


The suppression of noise is a subject of great interest 
because of its technological importance’ *. 

Elastic properties of the medium determine the 
mode of propagation and the speed of sound through 
it. In an un-ionized gas, the elastic forces are caused by 
the gas pressure which results from the collisions 
between the gas molecules. However, when a sound 
wave travels in a plasma, the electrical neutrality of a 
small element of plasma is disturbed owing to the 
relative displacement of the charged particles and 
= rau are called to restore the original state; 
; the elastic properties of plasma enabling it 
n acoustic wave. 


cts 


Ten eke 
‘ 


modes are possible, when an 
ough a plasma*. One of these 


ee = aM 
Saaz wave, in which the 


Be ee Oe a 


—— ms Serge = tol died 


where n is the frequency of acoustic signal travelling 
with a speed C, in air and C, in plasma and / the length 
of the plasma slab. x =[r,,+(1/r,.)], where r,5(p2C,/ 
p,C,) is the relative acoustic characteristic impedance 
of medium 2 (plasma) with respect to medium | (air) 
and Pp, and p, are the densities of the media | and 2 
respectively. If the two mass densities are equal, r, > 
=C,/C, and is equal to x provided r, » is large. It is, 
therefore, apparent from Eq.(1) that the ratio A3/A, 
depends upon n, /, the speed of IAW in the plasma and 
the speed of sound in air. 

If oblique incidence of the acoustic wave is 
considered and if C, is less than C, which is usual in the 
present case’ 7 (C,~5-10C,), then for some critical 
angle of incidence 0., the refracted ray makes an angle 
of 90° with the normal! to the interface: the critical angle 
(8.) is given by 


sn SC, /C; ~0.2-0.1 


Therefore 0. is of the order of a few nae If the 
angle of incidence is equal to or greater than a few 
degrees, no acoustic energy is transmitted to the second . 
medium and it is reflected back entirely 1 in the first 
medium. Thus, a very narrow cone of acousti riage es 
of significance if the transmission from air top plasma i 7 
considered. _ ee 

Numerical evaluation of si 
carried out by a computer 
C,. Some of the results are p re: 
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Fig. 2—Square root of sound power trans 


mission coefficient versus 


relative acoustic characteristic impedance of plasma (n= 5k Hz; /! = 1- 
10cm; C,;=5C,) 


Fig. 3—Square root of sound power trans 


relative acoustic characteristic impedance of 
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potentially powerful addition ¢ 


effecting analysis of a plasma and 
as sonic probe technique. 
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Fig. 4—Noise spectrum of a turbo jet engine: Pressure level versus 
frequency 


Noise has an essentially continuous frequency 
spectrum. The manner in which the intensity varies 
with frequency may be represented by a graph showing 
the sound spectrum of the noise. The abscissas of the 
graph are the frequencies and the ordinates give the 
pressure spectrum level of the sound at each frequency. 
The sound spectrum of a turbo jet engine of 150dB 
overall pressure level has been shown in Fig. 4. The 
spectrum extends from a frequency of 40 Hz to 10kHz. 
Thus turbo jet engines are quite rich in high frequency 
components which need to be suppressed. 

The above analysis shows that if such a noise source 
is surrounded by a plasma barrier of few centimetres 
thickness, the latter can effectively shield most of the 
noise spectrum from spreading further. 

The authors are grateful to Dr C P Gupta, Scientist 
In-charge, Computer Section, NGRI, Hyderabad, for 
carrying out the numerical analysis. They are also 
thankful to Department of Science and Technology, 
New Delhi, and to the University Grants Commission, 
New Delhi, for liberal grants, and to the authorities of 
Ravishankar University, Raipur, for the facilities 
provided. 
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A study has been made of resistivity (p), activation energy (AE), 
temperature coefficient of resistance (TCR) and mean free-path (/9) of 
charge carriers in air-oxidized antimony films in the temperature 
range 50-240°C. The electrical parameters of these films differ 
considerably from those of vacuum-evaporated films of antimony 
oxide. 

Dissociation and formation of new species has been 

observed during vacuum deposition of several oxides 

and chalcogenides' *. To avoid such dissociation, 

_ antimony oxide films were formed by oxidizing flash 

evaporated antimony films in air. Electrical properties 

— Of waGmony oxide in bulk have been reported 

; earlier’. Detailed studies on semiconducting, 

structural and optical properties of vacuum-deposited 

antimony oxide films*~ '? have already been made. In 

sonal eee note, an account is given of our 

‘es sstigation on semiconducting properties of 

1y oxide films formed by oxidation | of flash 
5 eons films. - 

a ASRS Is of the preparation of antimony oxide 

“films ‘os imila 10 those for eee oxide films. 

itimony was specpure. The domed of 


a aucun 
CG ie 


where ¢ is the electronic charge, E, the semiconductor 
energy gap, r the average linear dimension of the island, 
and K the dielectric constant of the material. Even 
though the films studied by us were much thicker and 
more or less continuous, the variation of AE has the 
same trend as that of very thin films. 

Temperature coefficient of resistance (TCR) was 
negative for thin as well as thick films. For each film, 
TCR decreased with increase of temperature. Higher 
value of negative TCR at lower temperature is also 
explained on the basis of the island structure theory. 
The expression for TCR is given by: 

dy babes 4 |. / 
Negative TCR = ‘Vo ) ie a) Tr? 


ns. a eae 


where C = [2e/er + AE,]/2K, d is the average distance 
between islands, m the electronic mass, € the dielectric 
constant of the substrate, @ the potential barrier 
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‘ 
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between islands and h the Planck’s constant. In the 
present investigation, the first term in Eg. (3) 
predominates. 

Mean free-path of charge carriers at different 
temperatures was obtained graphically for Sb2O3 films 
from pd versus d curves as done previously by 
Goswami et al.!’ Fig. 2 shows pd versus d curves for 
different temperatures and Table 2 shows the variation 
of mean free-path of charge carriers with temperature 
for a hypothetical massive material of Sb,O3 having 
the same structure and roughly the same number of 
lattice defects as those of the films. Thick films exhibit 
bulk or nearly bulklike properties. In a bulk 
semiconductor, [> increases with decrease of 
temperature due to lattice (intrinsic) impurities!*. In 
our present investigation, /p increases with decrease of 
temperature due to the presence of traces of SbO as 
intrinsic impurity in thick Sb,O; films which tend to 
behave like bulk semiconductor. 


1x10! 


10" 


dx 10°,¢m 
Fig. 2—Plot of pd versus d for different thicknesses 


= enmeeeremereee 
Table 2—Variation of Mean Free-path (lo) of Charge 
Carriers with Temperature 


Temperature (Io) 
"e 10°* cm 
100 1.224 
150 0.7823 
200 0.7442 
a 
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Antimony oxide films are expected to obey an 
analogue of Mattiessen’s rule. Hence the effective 
mobility of the film (Hg) Is obtainable from the sum of 
the inverses of the various mobilities, so that 


ee . (4) 
He Hy Hy = Hs 


where sl; is mobility due to lattice impurities, 4, that 
due to imperfections and Hs that due to surface 
scattering. 

Thus we conclude that films of antimony oxide are 
semiconducting in nature irrespective of their 
thickness. This result is only in qualitative agreement 
with the prediction of Neugebauer’s theory since 
during the gradual oxidation of antimony films, SbO 
was formed as an intermediate product which acted as 
lattice impurities in Sb,O, films. 

The authors wish to thank Principal R V Khairnar 
for his constant encouragement and keen interest in 
this study. 
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The effect of magnetic field (B) on drain current of FET BFW-61 
has been studied. A decrease in drain current on the application of 
magnetic field has been observed in all the three cases studied, 
namely B parallel to current flow (Bll), B perpendicular to the current 
flow and along the width of the crystal (B:\W) and B perpendicular to 
the current flow and along the thickness of the crystal (BiT). The 
change in drain current on the application of B has been found to be 
in the order when BT > B|| = BiW. The results have been explained as 
due to the combined effects of Hall effect and magnetoresistance. 


The characteristics, of p-n junction devices are 
susceptible to both temperature! ° and magnetic 
field®~'*. Field effect transistor is being used 
extensively since it combines the desirable features of 
the thermionic valves and those of conventional 
transistors and hence is a superior substitute for both 
in many circuit applications. In view of this, it is 
essential to know the degree of susceptibility of the 
characteristics of these devices to external influences 
like temperature and magnetic field for proper design 
of the electronic circuits. 

Studies on the effect of temperature on the transfer 
characteristics of a large number of samples of BF W-10 
have been carried out by Jamwal and Peer*. Studies on 
the effect of magnetic field on the drain current and 
transconductance of JFETs, BFW-10 and BFW-1I 
have been carried out by Jain and Agarwal’. In the 
present work, studies on the effect of magnetic field on 
the drain current of a large number of samples of BF W- 
61 have been carried out with a view to see as to what 
extent the characteristics of equivalent transistors are 
susceptible to the influence of magnetic field. The 
results are presented in this note. . 

The effect of magnetic field on the drain current of 
BFW-61 has been studied using the cirguit similar to 
- that employed by Jain and Agarwal” and by housing 
- the FET between the polepieces of an electromagnet 
(field strength was determined using Hall probe) and 


. ving the field along the three different directions 
sei current flow (Bll), (ii) B 
t flow and along the width 


crystal (BiW) and (iii) B perpendicular to the 


current flow and along the thickness of the crystal 
(BiT). These directions were determined by decapsulat- 
ing the FET and observing it under a microscope. 
The results on the measurement of decrease in drain 
current (AJ) with magnetic field applied in different 
directions for various gate-to-source voltages (Vs) are 
presented in Table |. Typical results when Vp>=15 V 
and Vg;=0V are shown graphically in Fig. 1. 
An examination of the data presented in Table 1 and 
Fig. | clearly shows that the decrease in drain current 
with field when applied in the BiT direction is more 
than the decrease found when the field is applied along 
the other two directions, viz. BiW and Bll. But the 
magnitude of change in drain current within the limits 
of experimental error is same in BiW and Bil directions. 
However, Jain and Agarwal® observed a greater 
magnitude of decrease in Ip on the application of 
magnetic field in BiW than in Bll directions for JFETs, 
BFW-10 and BFW-11 for a field strength of 19.5 kOe. 
FETs being unipolar devices and since magnetoresis- 
tance effect is a majority charge carrier effect, there is 
every possibility for magnetoresistance to appear in 
FETs. It is this magnetoresistance effect that is 
responsible for the decrease in the drain current when 


Table 1—Variation of Alp (in A) with B (in kOe) at 
Different Gate-to-Source Voltages for BF W-61 


(Vps=15 V] 
B Vos=0 V Ves=—0.5 V Ves=—1.0V 
(K Oe) a 
B\| Biw BiT Bil BLW Biw Bil Biw BT 


4.8 09 O08 10 05 G35 G6 Gb ip wa 
10.3 26° 25 40 “20° 20° 26 14 See 
15.4 ao 38 #49 #30 28 45- 49. 26225 
16.3 R22 4) 3S. 32 Bh, SY 2) eee 
17.3 50 $0 72 40: 39 GO 30. 2 ae 


B,k Oe 


Fig. 1—Variation of change in drain current (Alp) 
with magnetic field (B) for JFET BF W-61 
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the field is applied parallel to the current flow. In the 
other two directions, current flows perpendicular to 
the direction of magnetic field and therefore Hall effect 
is associated with the device. Hence application of 
magnetic field leads to the development of Hall voltage 
which makes the reverse biased gate junction to be 
more negative. Increased reverse bias increases the 
channel resistance leading to a decrease in drain 
current. However, the magnitude of decrease in drain 
current is found to be different in different directions. 
This behaviour suggests that the effect of magnetic field 
on the drain current in the FETs may be considered to 
be due to the combined effects of Hall effect and 
magnetoresistance. 

The studies carried out on a few other samples of 
BFW-61 FETs also indicate similar type of behaviour 
presented in this note. Even though BF W-10, 11 and 61 
are considered as equivalent in circuit applications the 
effects of external influences like temperature and 


“5 edi 


magnetic field are not alike on the characteristics of the 
FETs. This aspect is to be borne in mind in the 
application of these devices. 2 
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Electrical conductivity of naphthalene is measured at different 
voltages before and after its magneto-electret formation. It is 
observed that electrical conduction decreases by 5-20 %, depending 
on the forming parameters of magneto-electret. The variation of 
percentage decrease in electrical conduction with forming 
temperature closely follows the variation of percentage increase in 
dielectric constant with the forming temperature. An attempt is made 
to correlate the two phenomena. 


The observations on naphthalene samples showed’ 
that there is an increase (up to 5%) in its dielectric 
constant due to magneto-electret formation. On the 
other hand, measurement of current, which is being 
reported in this note, shows that electrical conduction 
in this material decreases (5 to 20%) on magneto- 
electret formation. The decrease in_ electrical 
conduction is an evidence of the loss of charge carriers, 
whereas an increase in dielectric constant indicates 
that the charge carriers get associated in such a way 
that the new configuration (which may be psuedo- 
dipoles or traps) has a higher polarizability. Since the 
two phenomena involve charge carriers in the bulk of 
the material, it is worthwhile to seek a correlation 
between them. The present note reports the results of 
measurements on electrical conduction of naphthalene 
magneto-electret and tries to correlate the obser- 
vations with those of dielectric constant reported 
earlier’. . . 
Magneto-electret samples of naphthalene (diam 1.25 
cm and thickness 1.25 mm) were prepared following 
the usual method? at different temperatures (30, 55, 60, 
65 and 70°C) using various magnetic field strengths 
(3.0-17.0kOe). The total duration of application of the 
| ic field was 1.5 hr (1 hr at constant elevated 
ture and 0.5 hr during cooling). 
ee caveat measurement was done at room 
temperature with the sample mounted in a holder 
similar to that described elsewhere*. ECIL electro- 
meter (EA 814) was used for the measurement of 
current at different voltages (in equal steps up to 500 V) 
slied to the sample before and after magnetic- 


electret formation. 


Fig. | shows variation of current measured at 
different voltages for magneto-electret (ME) samples 
prepared at the forming temperature of 60°C, using 
different magnetic field strengths. It is seen that current 
after thermo-magnetic treatment always has a value 
lower than that for the virgin sample. The decrease in 
electrical conduction ranges from 5 to 20% in most of 
the cases depending on the forming temperature and 
magnetic field strength. The difference between the 
values of current for the virgin and ME-samples 
increases systematically with the increase in the 
measuring voltage, though in some cases, especially in 
samples formed at low magnetic field, this difference 
seems to be constant. 

A systematic variation is observed on plotting the 
percentage decrease in current against forming 
temperature for various forming magnetic fields (Fig. 
2a). All the curves exhibit a maximum at the forming 
temperature of 60°C. Incidentally these curves have a 
striking resemblance to the variation of the percentage 
increase in the dielectric constant with forming 
temperature (Fig. 2b) which also exhibits a 
maximum at 60°C. However, the resemblance is only 
qualitative. The increase in dielectric constant ranges 
from 2 to 5% whereas the decrease in the current is 
from 5 to 20% and in some cases it is more than 20%. 
Yet, the qualitative resemblance between the two types 
of curves prompts one to seek for a correlation between 
the decrease in current and increase in dielectric 
constant caused in electret formation. 

The I-V characteristic of samples of naphthalene 
indicates that it perhaps follows the log/-V '/? law. In 


CURRENT, 10 ''a 


200 300 400 500 * 100 200 300 400500 
VOLTAGE ,V 

Fig. 1—I-V_ characteristics of magneto-electret (ME) samples 

prepared at the forming temperature (T,) 60°C using different 

magnetic field strengths (Broken curve: Virgin sample; Continuous 


curve: ME sample] 
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Fig. 2—Comparison of percentage variations of conduction current 

and dielectric constant due to electret formation of naphthalene 

samples with forming temperature [(a), conduction current decrease; 
(b), dielectric constant increase] 
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_ Fig. 3—Plot of log! versus V'? for magneto-electret sample 


prepared at forming temperature (7;) of 60°C using magneticeield 
strength (H,) of 3.0 kOe 


Fig. 3, log! 1s plotted against V'/?. At lower voltage 
ranges (up to 150 V), this curve is a straight line, 
whereas at higher voltages the curve transforms to 
another straight line of lesser slope. It seems that at 
lower voltages, the Poole-Frenkel process is 
responsible for electrical conduction, and this 
‘gradually changes to the Schottky process at higher 


In the Schottky process, the current is given by® 


I= AT? exp(—o/kT )exp(B,E"?/kT) - () 
where E is the applied electric field and 
B,=(e?/4ne96,)"? = a/e}!? .. (2) 


where a=(e*/4é9)'' is a constant and other symbols 
have their usual meaning. Eq. (1) shows that at a fixed 


™ 


temperature, log! o E", as has been observed. 


Further, if there happens to be any change in the 
relative dielectric permittivity (e,), it will have its effect 
on the current J. | = 
By differentiating Eq. (1) and using Eq. (2) it is easy to 
show that 
— dI/1 =(a/2kT)(E/e,)*? dé,/€, ae 
For a constant temperature and field, this relation 
reduces to 


—di/I x dé,/é, - 4) 


which shows that a fractional change in current is 
proportional to a practical change in the dielectric 
constant. The negative sign shows that as the dielectric 
constant increases, the current will decrease. This is 
what has been observed. 

It can be seen that Eq. (4) is applicable also for the 
case where conduction is governed by the Poole- 
Frenkel process. The Poole-Frenkel current is given 
by? 


1=04("F) | (Fie-)-* Jom Fer) ae (5) 


where B =(e?/mé9¢,)!/? = b/e}/? ... (6) 


where b is a constant. 


The measuring field E is of the order of 10°V/cm and 
the value of B/kT at room temperature is of the order 
of ten’. Therefore, 1 in the expression [(BE'/?/kT)—1] 
is neglected in comparison with BE'!?/kT. 

With this approximation, one can get on 
differentiation of Eq. (5), and using Eq. (6), 


—dI/I = 3[b/kT)(E/e,)'? — 1]de,/e, r 
d 


 < 


which at a fixed field and temperature again reducesto 
Eq. (4). a = 
Naphthalene magneto-electrets develop positive — 
isocharge*. This means that in contact with the metal cs 
(tin) this material loses some electrons to the meta 
elevated temperatures. This loss of condu 
electrons reduces its conductivity on subsequen 
measurement at room temperature. A fraction o 
electrons, while migrating to the metal at an 
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Electroluminescence in Cu-Phthalocyanine 
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Voltage and frequency dependence of electroluminescent 
brightness have been studied for Cu-phthalocyanine at room 
temperature. The brightness varies according to the relation B 
=B,exp(—C/V"”) in the low-voltage region and according to 
B=Byexp(—C/V) in the high-voltage region. 


The study of electroluminescence (EL) in organic 
materials has gained attention only in recent years’. 
Although EL in these materials is similar to those in 
inorganic materials, the fields necessary are much 
higher, the mean acceleration paths are much smaller 
and correspond to molecular size in such materials?. 
Helberger* observed the luminescence of magnesium 
and zinc phthalocyanine in boiling tetraline containing 
small amounts of benzyl peroxide, hydrogen peroxide 


or air. The present work is a study of dependence of the | 


EL brightness on the voltage (V) and the frequency of 
me spied electric field in Cu-phthalocyanine. 
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basis that initially the number of particles in which EL 
takes place is small, but on increasing the voltage, more 
and more active regions are exposed to voltage 
gradient above the threshold level. It is evident from 
the linearity of the plot of the mean EL brightness (B) 
with 1/V'/? that the relation between B and V is of the 
form: 

B= Byexp(— C/V''?),where Bo and C are constants. 
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NOTES 


case also, a slight deviation from linearity is found at 

higher voltages which may be because of variation of 

the type exp(—C/V). The increase of EL brightness 

due to increasing frequency above ~ 100 Hz is obvious 

because of the more rapid emptying and refilling of 

luminescent centres’. But when the time period of the 

applied ac field becomes comparable with the lifetime 

of the excited electron, the electron does not have 

enough time to emerge out from the traps; hence 

saturation can be expected in the high frequency range 

Curve Voltage,V which has been experimentally observed. The 

bo (PO decreasing nature of brightness at about 100 Hz can be 

understood on the basis of carrier injection and 
polarization effects®. 
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Efficiency of binderless electroluminescence (EL) as a function of 
voltage of the applied field is investigated. It is found that efficiency 
first increases, passes through a maximum and then decreases. A new 
interpretation of binderless EL in terms of Schottky-emission type of 
transport process is proposed. 


Electroluminescence (EL) has been investigated in two 
ways':(i) by studying the luminescence characteristics, 
and (ii) by electrical measurements. The first approach 
has been attempted by many workers and it has served 
as the basis for the mechanism of EL developed by 
Georgobiani and Fok’. The second approach is 
concerned with the measurement of absorption of 
energy of external excitation by the phosphor crystal; 
the efficiency of the crystal is defined as the ratio of light 
Output to the power absorbed by it. 

In intrinsic EL the phosphor is embedded in some 
dielectric and then pressed between iwo conducting 
plates one of which should be transparent. However, 
Bhushan’ has shown that if the electroluminor 
(ZnS,ZnO):Cu,Cl is treated with salt solution of 
CuSO, or Na,SiO;, EL can be observed just by 
pressing between two conducting plates (i.e. without 
the use of a binder or dielectric). Such a process is 
known as binderless EL. In the present note, efficiency 
of binderless EL has been investigated as a function of 
applied voltage for (ZnS,ZnO): Cu.C! electroluminor, 
from the ratio of light Output to the power absorbed. A 
new interpretation of binderless EL in terms of 
Schottky-emission type of transport process has been 
proposed. 

The electroluminor was prepared by firing a mixture 
of 10g ZnS (luminescent grade), 2g ZnO (luminescent 
grade), 0.093g CuSO,, 0.058g NH,CI and 0.2 
sulphur (all the three AR grade) at 1000°C for 2.S5hr. An 
extra amount of sulphur was poured on the upper part 
of the vessel which was then completely closed. The 
fired content was washed with 20°, acetic acid (to 
remove excess ZnO) and distilled water, dried and 
washed again with 4% NaCN (to removeexcess Cu) 
and dried. The electroluminor thus Prepared was 
treated with salt solution of Na,SiO, or CuSO, and 
then dried at about 120°C. 
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The EL cell was prepared by pressing the treated 
electroluminor between two conducting glass plates 
separated by 300 um thick mica sheet having a circular 
hole of area 2 sq cm for the electroluminor. For EL, 
excitation source was a transformer (50 Hz, 3kV). An 
RCA IP 21 photomultiplier (PM) tube, operated by a 
regulated power supply at 1000 V was used to measure 
the light output. The light output in the form of current 
was measured by a polyflex galvanometer (10° ® A/cm). 
The current through the EL cell was measured by a 
Philips GM 6009 VTVM by noting the voltage 
developed across a resistance of 1 MQ. 

Fig. 1 shows the variation of log of EL brightness (B) 
as a function of the inverse square root of voltage (V) of 
the applied ac field at 50 Hz for the electroluminors 
treated with Na,SiO; and CuSQ,. It is seen that except 
at lower voltages, the nature is well described by the 
relation B=Byexp(—b/V'!?), where By and b are 
parameters which depend on temperature, frequency 
of the alternating voltage, the phosphor and the details 
of construction of the EL cell. Since in the present case 
the studies have been done at constant frequency and 
temperature and the same cell has been used for all the 
studies, these parameters can be taken as constant. 

Fig. 2 shows the variation of current (I) through the 
EL cell measured by the VTVM under the excitation of 
ac field (50 Hz) as a function of applied voltage. From 
the straight line nature it is clear that [oc exp(aV"/2), 
where a is some constant. The efficiency (p) of EL has 
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Fig. 1—Variation of log B with applied 
(v8) voltage 
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been determined from the ratio of light output 
(expressed in arb. units) to the power input (expressed 


in watts) and plotted as a function of voltage (Fig. 3). band which then collide with the electrons in the 
With the increase in voltage the efficiency first increases luminescent centres which get excited and after return 
and then, after passing through a maximum, decreases. to the centres give rise to luminescent emission. The 
The variation of EL brightness according to the original electrons with their small amount of energy left 
relation B=Byexp(—b/V'?) indicates that the moveinthe crystal and again during the next half of the 
mechanism of excitation is of acceleration-collision ac cycle get ready for acceleration. Thus the higher the 
type, and a barrier of Mott-Schottky type formed in voltage of the applied field the higher is the light 
which the maximum field strength varies as V'/? Output; but for high voltages the power input increases 
faster than the light output* and hence the efficiency is 
maximum at some intermediate value of voltage. 
According to Yu L’vova', the value of EL efficiency 
is governed by the relative roles of (i) number of 
electrons participating in EL, and (ii) the conditions 
which exist during their acceleration; factor (i) is the 
dominant factor. The absorption of energy is governed 


(Ref. 4). According to this mechanism, electrons from 
donors or traps are first accelerated in the conduction 


and from the deeper traps. Thus by considering a 
kinetic of EL which consists of transport equations for 
a crystal phosphor with three levels, namely, 
luminescent centres, shallow electron traps and deeper 
donors, Yu L’vova has shown that EL efficiency (p) can 
be represented by the formula: 


B 1 
AV) = = V> exp(b, V+b/V') a 


ae vi/2 (vot) assumed that liberation of electrons from the local 
a Fig. 2—Variation of log! with applied voltage (V") keno pei Se Suneel Sanaa: ee 
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by the electrons liberated from the luminescent centres 


es ee : =; Where W represents the average power absorbed bya 
es 5 © 2 ee phosphor. In deriving this formula, it has been 
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Table 1— Values of Different Constants 


Treated «By x 10 es ax 102 
with Arb. units co LA A/V'!2 
CuSO, 6.4 124.368 5.4 1.302 


Na,SiO, 20 --. 158.304 3.3 7.992 


ecuSO%4 
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Log(W/v) 
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Fig. 6—Variation of log(W/V) with be Hage wage) 


it is assumed that the eouchectionn process is of 
Schottky-emission type, the current is ¢ 


1=T? [exp(—C+DV!)/kT] 


where C and D are constants. At ert tem 
this expression can be ak ti a 


I= Aexp(av’) “th sok fe 


where A and a are const ant: Tt us ti = ? 
be written as: a ee im 


< 
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NOTES 


os . " : > 
hird process is more effective in the binderless EL of 


(ZnS,ZnO): Cu,Cl. Eq. (5) has been found to be suitable 
yr CaS system also°. 

eather, the log(W/V) versus V curve (Fig.6) 

Be sovitce from straightness at lower voltages which 

Jicates that, at these voltages, shallow levels in this 

__ electroluminor participate in energy absorption’. The 


source of carriers may thus be taken as shallow traps or 
osete along with process (iii) as mentioned above. 


The author is thankful to Prof. RK Thakur for 
valuable discussions. 
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The electronic spectra of complexes of praseodymium(III) and 
neodymium(III) with amino glycine, alanine and valine have been 
recorded in the 360-920nm region in water. The energies and 
intensities of various transitions calculated using the Judd-Ofelt 
relation are in good agreement with those obtained experimentally. 


The absorption spectra of complexed rare earth ions in 
the f +f transition region have been studied by various 
workers! ’. Other spectroscopic studies have yielded 
information regarding the structure and composition 
of the complexes*~'?. In this note, we report the 
electronic spectra of neodymium(III) and 
praseodymium(III) complexes with glycine, alanine 
and valine. 

Various intensity parameters such as the Judd-Ofelt 
parameter (7,), the bonding parameter (b'/2), the 
energy parameters (F,, F 4, F,) and the covalency of the 
complexes have been calculated and discussed. 


Experimental Details 

The praseodymium(III) and neodymium(III) 
complexes were prepared by the method reported in 
literature’*. The chemicals used were of AR grade and 
the metals were 99.9 % pure (Indian Rare Earths Ltd). 
The complexes were synthesized with glycine, alanine 
and valine in molar ratio 1:1 and their formulas are 


[M(H,NCH,COO)8H,O]Cl,, 
[M(NH 7CH 7CH 2COO).8H 2,0 JCI 2 and 
CH; NH, 
Tespectively, where M represents the metal. 
The absorption spectra of the complexes were 
recorded on Carl-Zeiss VSU-2 Spectrophotometer in 
the 360-920 nm region in water solvent. The energies 


recorded in terms of wavelength were converted into 
wavenumber. 
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The results obtained for different complexes of 
Nd(III) and Pr(III) are. collected in Table 1. The rms 
deviation found between experimental and theoretical 
values are small, supporting the validity of the relation 
used. The red shift of the energy levels from that of free 
ion in different solvents supports the complexation**. 
The hypersensitive transition in the Nd complexes 
(419). > *Gs/2) shows a red shift when the ligands are 
changed in the order: alanine < glycine < valine. 

While in case of the Pr complexes for the associated 
pseudo-hypersensitive transition (*H, > 3P,) we have 
no such trend due to assumption of 4/-wavefunctions 
to be hydrogenic. It is observed that there is not very 
much shifting in peak position when ligands are 
changed, which may be attributed to some 
participation of oxygen atom orbitals only with those 
of the metal at different distances. Thus, the amino acid 
coordinates through carboxylic oxygen only while 
nitrogen remains unaffected. 


Table 1—Observed and Calculated Values of Energy (in 
wavenumbers) of Various Transitions in Nd and Pr 
Complexes 


NdA 


Complex* 
Levels 


NdG NdV 


NOTES 


Table 2—Observed and Calculated Values of Oscillator 


: Ci : Table.3—Computed Values of Intensity Parameters T, (in 
Strength for Various Transitions in NdA*, NdG* and NdV* em *) for Nd** and Pr?* Complexes 
Complexes Complex T, x 10° 
Complex* NdA NdG NdV ss 
= : = T, T, T; Ts 
Levels Exp. Gk Exp. Cale. Exp... Cale, NdA 0.438 0.7355 0.6469 1.137 
x 10 x 107° x 107° NdG 0.524 0.4775 0.7163 0.6666 
2P. 0.639 0.658 0.177 0.426 0075 0.430 NdV 0.3995 0.4835 0.5505 0.8782 
a 0.876 0.198 0.163 0.179 0.315 0.152 PrA ~ 4.057 1.05 3.165 0.332 
*G.,, 0.765 1.567 0.174 1.287 0.144 1.138 PrG 2.63 1.2445 3.27 0.380 
*Gor 0.744 0389 0.556 0.336 0276 0.291 Prv - 6.769 0.637 2.671 0.2385 
"Gus, 3.387 3.364 2.567 2.749 2.481 2.461 Se eee ee - ——s 
7 
403 12731 12.485 11.576 12.011 9.740 10.010 Table 4—Computed Values of Bonding Parameter f (in cm~') 
Fo) 0.450 0472 0.208 0476 0.181 0.382 B''? and b'? for Nd?* and Pr?* Complexes 
‘F,. 3.932 4099 4.111 4344 3.154 3.413 ‘ 
*F,,  SATL 5368 4940 4941 4084 4.134 Cam giek = =f, Fe Fe B As 
"Fe. 2.163 2356 1.790 1.710 1.333 1.624 Free ion 331.159 50.714 5.154 = = 
rms Dev. +0.3742 +0.3730 +0.3680 NdA 333.58 S1.917 = 5.428 ~—- 1.0073 — 
= ine eas NdG 336.285 50.717 5.494 1.015 ae 
NdV 336.927 48.317 5.500 ‘1.0174 a 
Esp: Cole, Erp. .Cale .. Exp. Cale. ax ga ia hn 170 0.9668 0.1288 
3p, 10.54 10.54 11.046 11.046 8694 8694 po 31133 4298 470 09665 0.1290 
5p, 6.718 3.838 6.26 4.549 4.156 2.328 PrV 311.40 42.98 4.70 0.9668 0.1288 
3P, 0.944 3.772 2.790 4470 0.498 2.288 
'D, 2.924 2924 3.362 3.362 2846 2.846 
rms Dev. +2.018 + 1.198 + 1.279 


J 
= 


Note regarding abbreviations A, G and V the same as in Table 1. 


Intensity parameters—The experimental values of 
the oscillator strengths for the various transitions have 
been computed by using the relation: 


P=46x 10° x¢,, x Av, 


where é,, and Av,,, are the molecular extinction 
coefficient and bandwidth for transition respectively. 

Using the Judd-Ofelt relation in terms of 7, 
parameters, we have 


Poa =[U Pv, T> + (UV ¥.T, + [UO VV. T 


and matrix elements [ U“ ]? represented by Carnall et 
al.'*, we calculate the oscillator strength for the f+ f 
transition (Table 2). The rms deviation between the 
observed and calculated values is very small. The 
parameters T,, 7, and T, obtained for the complexes 
are collected in Table 3. The parameter 7, is sensitive 
to the environment but its negative value has no 
physical significance. 

It is observed that the ratio of 74/7. is almost 
constant for the same metal but is more when Nd is 
replaced by Pr. This supports the fact that the 
coordinating environment around the metal ions 
remains unaltered, as mentioned earlier. 

Bonding parameter—The parameter aa which is a 
measure of the partial covalency in a metal-ligand 
bond can be.calculated from the relation pb? = [(1 
— p)/2]"? and the results are given‘in Table 4. In case of 


ae 8 i 


Nd, 51/2 is negative indicating that the complexes are 
ionic. The nephelauxetic ratio B together with F,, F4, 
F., for the Nd and Pr complexes are also shown in 
Table 4. 

One of us (AKG) thanks the University Grants 
Commission (India) for providing a fellowship. We also 
thank Prof A N Nigam and Dr S P Tandon of the 
Department for helpful discussions. 
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As a part of the study of iron-bearing minerals in the eastern 
Himalayan foothills, the ferrous-ferric ratio for a number of 
weathered samples in agricultural and forest-covered lands in the 
altitude range 100-2700 m has been estimated by the Mossbauer 
method. The results have been compared with similar such results 
for other countries with a view to indicating the trend of behaviour 
of iron in bulk and separated samples. The importance of these 
results in the field of geology and agriculture has been pointed out. 


A study has been undertaken to investigate the iron 
content in rocks and minerals in the eastern 
Himalayan foothills. As a part of this eco- 
development programme, soil and rock samples from 
a number of localities have been analyzed for their 
iron content using Mossbauer spectroscopic method. 
The Mossbauer spectroscopic procedure enables the 
analysis to be made in a non-destructive way and is 
able to distinguish between the two important 
oxidation states of iron (Fe? * and Fe3*) that occur in 
the structure of soil and rock under natural 
weathering. Accurate estimates of Fe?*/Fe3+ ratio 
are very important in the field of geology and 
agriculture. In geology, this ratio can be taken as 
indicative of the temperature and pressure of 
formation of minerals which constitute the various 
rocks. In agriculture, this ratio can be taken as a 
measure of the degree of iron deficiency and hence an 
index of applicability of fertilizers. The determination 
of this ratio from chemical analysis is subject to 
uncertainty due to difficulty in dissolving rock 
minerals and _ resulting Possible oxidation. 
Consequently, the estimates from chemical method 
do not reflect the actual chemical States of iron in the 
minerals. The Mossbauer method, however, provides 
a rapid and accurate method for determining this 
ratio. 

The weathered soil samples of A-horizon were 
collected in autumn of 1979 from localities in the 
range 22°N, 87°E to 27°N, 89°F at altitude levels in 
the range 100-2700m. Each sample under field 
condition was first disaggregated and the fine 
fractions of sizes less than 300 um were separated by 
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sieving. An electromechanical transducer operated in 
the constant acceleration mode coupled to Nuclear 
Data 1100 Multichannel Analyzer was used with a 
14.4 keV gamma-ray source (°’Co in palladium 
matrix and of 3.8 mCi strength). In all measurements, 
Mossbauer spectra of samples were taken at room 
temperature. Fe?* and Fe** states were assigned to 
doublet components in the spectra using the 
procedure of other workers’. In the spectrum of 
CTRI sample (Fig. 1), the doublet AA’ is assigned to 
Fe** and the doublet CC’ to Fe**. Fig. 2 exhibits 
Mossbauer spectrum of BDTS sample showing only 
the doublet assigned to Fe**. The sample CTRI was 
obtained from a tobacco-field where the annual 
rainfall is about 400cm. The sample BDTS was 
obtained from field under tea plantation where the 
annual rainfall is about 400 cm. 

The Fe’*/Fe** ratio in a sample has been 
determined from the area ratio in the following way. 
If a thin absorber sample contains iron atoms in two 
different chemical states which give well-resolved 
Mossbauer spectra of area A, and A,, then 4 ,/A2 
=n,/n,, where n, and n, are the total number of iron 
atoms of the first and second types respectively. This 
is true when the full width at half maximum of the 
Mossbauer lines corresponding to two types of iron 
and their recoil-free fractions are equal in a very thin 


4s ~ : 
VELOCITY (mm/sec ) 


Fig. 1—Méssbauer spectrum of CTRI 
(Points represent experimental obser- 
vations, the solid line gives the curve 
obtained by fitting the experimental 
points and the dashed curve gives the total 
envelope of the fitted curves) 


NOTES 


Table 1—Fe** /Fe?* Ratios for Different Samples 
Sample Country Fe**/Fe?* Reference 
No. 
| 907 India 0.9345 Present work 
: 908 do 1.3157 do 
; 909 do 1.2610 do 
| 910 do 2.1413 do 
911 Zaire 1.6393 Bruyneel er a/.5 
1027 do 0.8547 do 
2220 do 1.3333 do 
2320 do 1.6666 do 
2514 do 1.4925 do 
912 Israel 2.2222 Rozenson and 
Heller-Kallai? 
age ee ee 
-2 ) 2 Table 2—MéOssbauer Data for Different Samples 
. VELOCITY ( ma/sec ) : ; 
Fig. 2—Méssbauer spectrum of BDTS (Mossbauer spectra did not show ferrous doublets) 
. (Points, solid and dashed curves mean the Sample Country Fe3+ Fe?*/ Reference 
! same as in Fig. 1) No. Fe3+ 
IS* = QSt 
S; India 0.57 0.71 0.0 Present 
2 w (NDIA) wok 
S-ESLInGeN et af SENTONITE : Si do 052: 0:79 = 06-2 
N B, do 0.53 0.72 0.0 do 
ee re eer iwak BDTS do 040° = O75.2- aay. hee 
- a cy do 0:57 40:65 0.0 do 
M, do 0.64 0.56 0.0 do 
M, do 0.61 0.47 0.0 do 
aS donk OR.| Dat. 0.0 do 
: LT-42 Montana 0.36 — 0.69 0.0 Eslinger 
| ~ etal.* 
LT-139A — do -0,35 0:59 — 00 do 
: _ = —. LTE-1 do 0.37. 41.09 0.0 do 
: SS 5 <a Lt-1 do 0.35 0.60 0.0 do 
iS 


*IS, Isomer shift (in mm/sec) with respéct to sodium siete, 


_ §-BULK SAMPLE OF WEST BENGAL(INDIA) FQS, Quadrupole splitting (in mmysec) 


@ —-ROZENSON & HELLER -KALLAI?, 
- MONTMORI-LLONITE (ISRAEL) 


; " ; -s-eRurNedt et al? BULK SAMPLE (ZAIRE) 


_ together with the results in —— = 
samples from Israel indicate a trend of a Bt 
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A data processor interface for the indigenous hardwired 4K MCA 
has been developed. The system consists of Micro-2200 
programmable calculator with its peripherals, micro data interface 
and a floppy disk drive. The data interface is designed for the transfer 


of data from MCA to Micro and reverse transfer. The hardware of - 


the GET DATA and INV GET DATA functions is briefly described 
in the present work. The necessary software for programmed transfer 
of data and for algorithms on statistical analysis of spectra has been 
developed for the system. The system can be conveniently used to 
analyze standard nuclear y-ray and X-ray spectra. 


Multichannel analyzers (MCA) have been extensively 

used in experimental investigations since early sixties. 

With the semiconductor technology, their memory 

capacity has increased to 4K-8K channels. The use of 

microprocessors in the MCA has made such features of 
| data acquisition and on-line data processing possible, 
S which would be impractical to implement in hardwired 
a. MCA. 

____ In this note, we briefly describe a data processing 
extension to the hardwired 4K channel MCA designed 
by the Computer Section, Bhabha Atomic Research 

Centre. The data processor consists of a programmable 
pak , sage Micro-2200, its peripherals; micro data 
: (MDI), floppy disk attachment, magnetic 
d | reader (MCR) and X-Y recorderf and a data 
face specially designed for the transfer of data from 
At to Micro-2200 and the reverse transfer. For the 

nuclear or X-ray spectroscopy experiments, 
1 processing can be accomplished while other 
are accumulating. Spectrum analysis 
s-in tegration, peak search, net area and 
culation, energy calibration, transfer of 
can be easily carried out on the system. 
disk attachment the mass storage of 


” of the complete system is shown in 
unicates _ a devices 


simultaneously interfaced to MDI, with interrupt 
facility. The hardware of the data interface connecting 
MCA to Micro consists of two parts (i) GET DATA, to 
read the region of interest from MCA to Micro (ii) INV 
GET DATA, to send data from Micro to MCA 
sequentially. 

GET DATA interface—Fig. 2 shows the block 
diagram of the interface. The transfer is done by 
pressing the DATA key on the keyboard of Micro. On 
the execution of this key function, Micro reads data up 
to a maximum of 9 digits in steps of two digits at a time. 
The five lower significant digits are used for the channel 
contents and the next three digits for the channel 
address. 

When the DATA key is pressed, Micro generates 
SEND, DEVICE STROBE and ACKNOWLEDGE 
signals. Region of interest (ROJ) signal from MCA is 
also used in Data Interface. A complete description of 
the signals available from MDI and their waveforms is 
given in the HCL manual on MDI and the circuit 
diagram shows the link of MDI and data interface. 

Since Micro accepts 2 BCD digits, i.e. 8 bits at a time, 
and the data to be fed from MCA to Micro is 32 bits (8 _ 
digits), multiplexer chips of type 74153 are used for this 
purpose (Fig. 2). By pressing DATA key, transfer of one 
channel number and its content is done. Micro can be 


programmed to give the transfer of a number of ROIs % 
aan 


registered on MCA to Micro. a 
INV GET DATA Interface—The control logic of the | a 
hardware of the interface is similar to the logic used in 
GET DATA function. To get data on MCA, MDI 
accepts two digits at a time from Micro and with the — 
signals TAKE, DEVICE STROBE, 
ACKNOWLEDGE and DEVICE NO, sends them to. 
the interface. The two-digits input from Micro is mee 
in latches to produce 32 bits output. With the help of | = 
latches 7475 and four demultiplexer chips of ee = 


DEVICE NO. 2 5 oo 
2 DEVICE NO.3 3 - 
a= 
1 - ro. 
DEVICE NO4 om 
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Fig. 2—Block diagram of the GET DATA Interface 


74155, a 20 bits output, containing the lower significant 
5 digits of Micro register representing the channel 
contents, is stored on MCA. The channel increment for 
reverse transfer is sequential. 

SOFTWARE 

A software package of a number of programmes has 
been developed on Micro-2200 and is stored on floppy. 
A brief description of these programmes and 
subroutines is as follows: 

Transfer routines: This is the main programme for 
the programmed transfer of data from MCA to Micro- 
2200 and vice versa, with the GET DATA and INV 
GET DATA keys on the keyboard of Micro. The 
programme also includes storage of the spectrum on 
the floppy disk. 

Spectrum smoothing: Smoothes the spectrum 
stored on floppy disk using either three-point or five- 
point smoothing functions, and stores it back on the 
disk. 

Peak search: Finds the peak start, peak channel 
and peak stop in a spectrum of MCA or on floppy. 


Peak analysis: The peak regions of interest ROI are 
selected on MCA using the thumbwheel switches on 
the control panel of MCA. 


The programme does peak integration, background 
calculation by linear interpolation between the 
extreme channels, peak net area, FWHM and % error 
in area determination for each ROI and prints the 
results. 

Other programmes are developed for specific use to 
analyze dispersive X-ray fluorescence spectra. 

A complete listing of the programmes and the 
detailed circuit diagrams for the interfaces are available 
from the authors on requiest. 


We would like to express our thanks to Shri J D 
Gupta, Sudhir Jain, C G Chatur, B P Dube and G § 
Lodha for their help in designing the interface. We also 
thank Dr N S Satya Murthy, Dr S S Kapoor, Dr V§ 
Ramamurthy and Shri Virendra Singh for their 
guidance and help on various aspects during the 
development of the system. 
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Uranium concentration in water collected from several parts of 
Assam, India are determined in clear and stirred conditions by the 
fission track method. In the clear condition U-concentrations are 
found to vary from (0.08 + .002) to (5.32 + .02) yg/l. It is observed 
that the U-concentration increases in the stirred condition. The 
results are compared with the global background. 


Uranium is the heaviest trace element occurring in 
nature in almost all minerals, rocks, water, etc. It is, 
therefore desirable to have a knowledge about 
concentration of uranium in natural samples like water 
from various sources. Such a study may be useful in 
throwing new light on the environmental physics of a 
region and for health physics as well. Uranium 
concentration which varies from ppb to ppm level can 
be measured by different methods, the latest being the 
charged particle track analysis first suggested by Price 
and Walker in 1963. The advent of the inexpensive and 
simple solid state track detectors’ makes possible high 
sensitivity measurement of mineral concentrations of 
fissionable impurities or of #-emitters. In our present 
experiment, uranium in water collected from different 
parts of the state of Assam, namely Gauhati, Barpeta 
and Mangaldai sub-division, was measured by using 
solid state nuclear track detectors. 

The water samples from different sources were 
collected in plastic pots which were first washed with 
boiling 25% HNO; acid. From each sample, 
subsamples were taken by micro-pipette with and 
without stirring the original samples. In order to 
maintain the amount of stirring for various samples at 
the same level an ultrasonic vibrator was used. 

A known volume: of water (0.05 ml) using 
micropipette, rinsed with distilled water and then test 
water, was placed on a clean circular disc of lexan 
detector. The water droplet was allowed to evaporate 
in still air placed inside an oven (at ~ 70°C). As shown 
in Fig. 1, this droplet would leave a thin layer of 
nonvolatile constituents including uranium. 
Ordinarily, water droplets would evaporate from the 
surface of the lexan detector in the manner shown in 
Fig. 1. 

The diameter of the droplet remained same during 
the time of evaporation. Then the lexan was covered 
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with another piece of lexan detector. The lexan 
detectors along with a standard glass of known 
uranium concentration serving as flux dosimeter were 
exposed to thermal neutrons at the CIRUS reactor, 
Trombay. 

After.irradiation, the detectors were separated and 
washed. The lexan detectors were etched with 6N, 
NaOH solution at 70°C for 14 min and then were dried 
up in air. Thereafter the entire surface area of the 
droplets on the detectors was scanned. In order to find 
the thermal neutron dose for these samples, the 
irradiated standard glass was freshly cut and etched in 
20% HF at 20°C for 15 sec and then tracks were 
counted under a high-power optical microscope. 

The general formula for uranium concentration (C,,) 
in units of weight/volume of liquid is given by” 


Nm 


Os VGN ooEg 


where N is the total number of tracks counted on the 
lexan; V, the volume of the liquid drop; m, the atomic 
weight of uranium; o, the reaction cross-section of 
38; ~, the total neutron dose and No, the Avogadro 
number. 


G is a geometrical factor that accounts for the 
number of charged particles produced per neutron 
capture and the fraction of those that register in the 
track detecting material. Here the value of G is 1 since 2 
fission fragments are produced per reaction and half of 
these reach each detector plate if the uranium deposit is 
thin. 

The factor E=1—(V,/V,)? is the etching efficiency 
where V/V, is the etching ratio defined as the linear 
rate of dissolution of the detector divided by the rate of 
dissolution along the damage track itself. In the present 
case, as the factor E is close to unity, the dependence of 
C,, on E may be neglected. 

The uranium contents of clear and stirred water 
samples are given in Table 1. The difference of the 
uranium content in the clear and the stirred conditions 
was deemed to have arisen due to the suspended 
particles which might have settled down under clear 
THIN OROPLET OF 


NONVOLATILE 
CONSTITUENTS 


O-WATER DROPLET 


EVAPORATING 
WATER DROPLET 


L—LEXAN DETECTOR 


Fig. 1—Diagram showing the stages of evaporation of water 
droplets from the surface of the lexan detector 
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Table 1—Suspended and Dissolved Uranium in Water 
Sample Clear sample Stirred sample Diffe- 
rencet 
No. of U-conc. No. of U-conc. yg/| 
tracks beg /I tracks pg/l 
M —I 13315 0.70+.006 15777* 0.83+.007 0.13 
-il 3445* 0.18+.003 3873 0.20+.003- 0.02 
—IIl 3027* 0.16+.003 3226 0.17 +.003 0.01 
—IV 3655 0.19 + .003 — — . 
—V 6789 0.36+.004 13330 0.70 + .006 0.34 
—VI 1562* 0.08 + .002 — — — 
B-—VII 101798* 5.32+.02 126999* 6.64+.02 [32 
— Vill 1519* 0.08+.002 1123 0.06+.002 —0.02 
—IX 15967 0.84+.007 40690 2.13+.01 1.29 
—X 16947 0.89+.007 40133 2.10+.01 1.21 
—X] 27785 1.45+.008 78234*  4.09+.02 2.64 
—XIl 40421 2.11+.01 _ — — 
G—XIll 2977 0.16 + .003 — — — 
— XIV 6668 0.35 + .004 — — — 
—XV 3271 0.17 + .003 a aS cor 
—XVI 2064 0.11 +.002 — = — 
—XVII 6704 0.35 + .004 —- — — 


* Associated with fission stars 
tSuspended uranium in particulate matter 


condition. The uranium contents are given with one 
standard deviation derived from the number of tracks 
counted. In all but one case (B-VIII), the uranium 
content in stirred condition is found to be higher than 
that in clear condition. This observation (B-VIII) is 
perhaps due to the fact that the sample contains fission 
stars (due to particle in suspension) only in the clear 
state. 

The measurement of uranium content of water by 
the methods other than this method is (generally) 
complicated due to the presence of variable mass of 
suspended uranium-bearing particles in water. But the 
fission track method described here makes the 
measurement simple and less expensive. 
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Applying fission track method, Fleischer and 
Lovett? determined the U-content of city water 
from the Pasadena, California, water supply and 
found it to be 9 pg/l, which was roughly three 
times greater than the U-concentration normally 
found in sea water®. Fleischer and Delany* deter- 
mined U-content in nautral water collected from 
different sources and found it to vary from (0.08 + 0.01) 
to (5.27+0.28) ug/l in the case of clear water samples 
and from (0.32 + 0.02) to (7.18 + 0.14) ug/l in the case of 
stirred samples. Chakarvarti et al.” determined the U- 
concentration of water collected from hot springs, river 
water, lake water, etc. They found that the U-content of 
hot spring water samples from Puga valley, Manikaran 
and Vashist kund varies from 1.4 to 7.4yg/l; 2.2 to 
2.8 ug/l and 2.7 to 4.2 ug/l respectively. In the case of 
river water and lake water, the U-concentration as 
obtained by Chakarvarti er al.° indicated a variation 
from 4.6 to 6.6 pg/L. 

The range of U-content as obtained in this 
experiment lies within the general range of U-content 
for background water observed by earlier workers. We 
observe nothing abnormal in the U-content of various 
water sources of the region. 

The financial assistance from the CSIR, New Delhi 
for this project is gratefully acknowledged by the 
authors. The authors wish to thank the Heads of the 
departments of Physics and Geology, Gauhati 
University for extending necessary laboratory 
facilities. 
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initials of authors, title of the paper (underlined), name of the 
conference, and where and when it was held (e.g. Herczeg P, 
Symmetry-violating kaon decays, paper presented to the 
International Conference on High Energy Physics and 
Nuclear Structure, Vancouver, Canada, 13-17 August 1979). 

Reference to a thesis should include the name of the — 
author, title of the thesis (underlined), university or a 


institution to which it was submitted and year of submission 
(e.g. Mehrotra S N, Many-bod ques and their 
applications _to interacting bosons, Ph D thesis, Ranchi 4 


University, 1976). ~ a 
Reference to a patent should include names of patentees, Son 
country of origin (underlined) and patent number, the _ 
organization to which the patent has been assigned n 
circular brackets), date of acceptance of the 
reference to an acting periodi availat 
Labes M M, US Pat. 4,066,567 (to ’ empie University), . 
January 1978; Chem. Abstr., 88 (No. 20) (1978), 138350 n] 
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